mew STICA 


Vol. 3 


1953 


No. 3 


THE PRACTICAL REPRESENTATION OF STANDING WAVES 
IN AN ACOUSTIC IMPEDANCE TUBE 
by W. K. R. LIPPERT 
Division of Building Research, C.S.I.R.O., Melbourne, Australia 


Summary 


It is shown that the envelope curves of the standing waves along the whole of an acoustic 
impedance tube with attenuation can be represented with the aid of diagrams. Suitable diagrams 
for easy application to impedance measurements are given and their usefulness is discussed in 


some specific cases. 


The position of the pressure minima are measured with different orifices of the probe tube and 


the resultant shifts of position are discussed. 


Special test measurements with a rigid terminal are proposed for specifying the tube per- 


formance. 


Sommaire 


On montre que les enveloppes des ondes stationnaires dans un tube absorbant 4 impédance 
acoustique peuvent étre représentées a l’aide de diagrammes. On donne des diagrammes utilisables 
dans les mesures d’impédance acoustique, et on démontre l’utilité de la méthode proposée par 


quelques exemples expérimentaux. 


Les positions des minima des ondes stationnaires ont été mesurées avec différents orifices du 
petit tube sonde; on discute les variations correspondantes de position. 
On propose de faire quelques expériences propres a caractériser la qualité d’un tube acoustique. 


Zusammenfassung 


Es wird gezeigt, dafs die Umhiillungskurven der stehenden Wellen in einem akustischen 
Impedanzrohr mit Dampfung langs des ganzen Rohres mit Hilfe von Diagrammen dargestellt 
werden kénnen. Geeignete Diagramme fiir den praktischen Gebrauch bei Impedanzmessungen 
werden angefertigt, und ihre Niitzlichkeit wird an Hand einiger Messungen diskutiert. 

Die Lage der Schalldruckminima wird mit verschiedenen Offnungen des Abtastrohres gemessen, 
und die sich ergebenden Verschiebungen werden diskutiert. 

Besondere Priifmessungen mit starrem AbschluB werden vorgeschlagen, um die Giite eines 


Impedanzrohres zu kennzeichnen. 


1. Introduction 


The transmission line method is the most 
widely used one for measuring the absorption 
coefficient or the specific acoustic impedance of 
a material. A smooth rigid-walled tube has a 
source of sound at one end and the sample under 
test at the other and the acoustic impedance can, 
as is well known, be found by measurements of 
the standing waves along the tube. 

The method requires the measurement of the 
pressure in the tube without disturbing the sound 
field, and several approaches to the solution of 
this problem are known [1]---[6]. That of TayLor 
{1], improved by Scorr [5], uses a long smallbore 
probe tube to measure the sound pressure in the 
main tube. The probe tube is moved along the 
main tube and is terminated by a microphone 
outside the main tube. 

This method is frequently used because of its 
simplicity and accuracy. Scorr has shown that 


it is not permissible to assume that the tube 
attenuation is negligible when precise measure- 
ments are required and he derived the relevant 
theory. The practical application of this theory 
is lengthy, particularly when the standing wave 
has been explored by recording sound pressure 
by a high speed level recorder. BERANEK [7] dis- 
cussed Scott’s theory and showed how to derive 
the specific acoustic impedance from measured 
data from part of the standing wave close to the 
sample. However, in the course of the present 
work it was found preferable to use the record 
along the whole of the tube, and suitable diagrams 
as presented below were found particularly use- 
ful. With these diagrams a sensitive method is 
available for detecting irregularities in the oper- 
ation of the equipment, extrapolation of minima 
to zero distance is easier and the whole provides 
an improvement in the specification of tube per- 
formance. 
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2. The envelope function 


Acoustic impedance tube theory for a tube with 
attenuation is given in text books [7] and it is 
necessary only to quote the results here. The 
ratio of the Mth pressure maximum to the Nth 
pressure minimum is 


(Fe aew (1) 


2 cosh? (x Du + v1) 


ae 2 (« Dm + y,) 


_, sinh? 2 («dn + y,) 


2 sinh? (« dy + oan + = 


and the distance dy of the Nth pressure minimum 
(counted from the sample) is 


dy 2N—-1 yy 
A 4 2a 


er. sinh2(ax+y,), (2) 
wherein « is the attenuation constant for the 
sound wave due to energy losses at the side walls 
of the tube and in the gas, 


f=/22 is the frequency, 
c is the wave velocity, 
A is the wavelength, 
k=-w/c=2a/A is the wave number, 
x is the distance from the sample 
surface, 
dy is the distance of the Nth mini- 
mum from the sample surface, 
is the distance of the Mth maxi- 
mum from the sample surface, 
P; is the pressure of the incident 
wave, 
P, is the pressure of the reflected 
wave and y= y, + j is given by 


P,| Pi ee ea 7 a Peg ate iv.) (3) 
The values of y, and y, are calculated by sub- 


stitution of experimental values of ( Prax /Pmin) m,n 
and dy/A in equations (1) and (2) and the acoustic 
impedance ratio is then 
vie R Xx 1 
j—=- h | Wo) & 
ae inp eer coth (y; + jp.) 

Vk (4) 


~ coth (yp, + j p2)» 
wherein Z_ is the acoustic impedance, 
R_ is the acoustic resistance, 
X is the acoustic reactance, 
and 0 is the density of the medium. 


Du 


The ratio of the maximum to the minimum 
pressure in the tube is commonly measured by 
moving the probe tube at constant speed and 
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recording the sound pressure with a high speed 
level recorder. A typical record is given in Fig. 1. 


The curves (5) 
P= al/2 cosh? («x + Be — a sinh? 2 (xx + yy) 
and (6) 


P=A 2 sinh? (ax + 1) + 375 D aehe 2 (ax + Y), 


ok 
where A is a constant depending upon the sound 
pressure level, will enclose the whole of the re- 
corded trace and hence are envelope curves. 


x 


Surface of sample atx=0 


Oe curve 


Envelope 


Minima curve 


—>|104B|=— 
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Fig. 1. A typical record of the sound pressure with a high- 
speed level recorder. 


These curves have been superimposed on Fig. 1. 
The curve which encloses the maxima (cf.eq.(5)) 
is practically a straight line since «x is usually 
small. The envelope curves (cf.eq.(5) and (6)) 
will be separated by an ordinate (in dB) which 
may be called the envelope value L (x) such that 


L(x) = (7) 


2 cosh? 8 4 ”) —~ ae sinh? 2 («x + y,) 


20 log 
2 sinh? es + ray + seat sinh? 2 (ax + y,) 


which can be readily transformed to 


L(x) = (8) 


1— cs inl 2 (ox +) 


20 log | coth (xx + y,) 

1+. — mans (xx + wy) 
It will be observed that the shape of the envelope 
curve for a given tube and frequency depends on 
the tube attenuation and the impedance of the 
sample. The value of y, (the phase change caused 
by the sample) does not affect the envelope; a 
change in y, affects both the value and position 
of the maxima and minima in such a way that 
the envelope curves are unaltered. 

Before discussing the envelope curves further 
it is desirable to distinguish two cases: case A 
including highly reflecting to normally absorbent 
samples, say y, 0.6, and case B with highly 
absorbent samples, say y, > 0.6. 
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3. The envelope function for highly reflecting 
to normally absorbent samples (Case A) 


When y, $0.6, i.e. when |P,|/|P:| 20.3, the - 


envelope function L(x) can conveniently be 
divided into two parts 
L(x) = L,(x) — Ly (x) (9) 
where L,(x) = 20 log coth (xx + y,), (10) 
2 

le ae cosh? («x + yy) 
L,(x) = 10 log ; oe (11) 

1 —~ sinh? (ax + y,) 

ja sinh? (a Wy 


L,(x) may be called the principal envelope 
function because it largely determines the en- 
velope; L,(x) is the correction function, small 
compared with L,(x). The relation between the 
various functions and envelope curves is shown 
diagrammatically in Fig. 2. 


Measured - 
envelope 


Fig. 2. Diagram of relation between the defined functions 
and envelope curves. 

The principal envelope function L, (x) has been 
plotted against «x in Fig. 3 for various values of 
L,(0) and for ease of application a set of scales of 
the distance x for different values of « is given. 
In Fig. 4 the correction function L,(«) is plotted 
against («x + y,) for various values of the para- 
meter «/k; it will be observed that the value of 
L(x) depends primarily on the value of «/k. 

The attenuation constant « may be calculated 
with sufficient accuracy for a smooth rigid-wall 
tube when the size and frequency are known 
[7]---[12]. Beranex [7] recommends that the 
values obtained from KrrcuHorr’s equation be 
increased about 8 per cent in magnitude when 
predicting the attenuation and this gives an at- 
tenuation constant « in nepers per cm 


U 


a=1.59-10-*- fh. =, (12) 


wherein fis the frequency in c/s, 
U 
S 


For normal frequencies and dimensions, the at- 
tenuation constant « is a very small quantity and 
consequently «/k is normally much smaller than 
unity. 


is the circumference of the tube in cm, 


is the area of the tube in cm?. 
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Fig. 3. The principal envelope function L,(«) plotted 
against «x for various values of L,(0) and a set of 
scales for different values of «. 
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Fig. 4. The correction function L,(x) plotted against 
ax +y, for various values of «/k. 


If now it be required to determine y, corres- 
ponding to a particular curve recorded in the 
impedance tube, the values of L,(x) are found 
from the value of «/k and from Fig. 41 and the 
curve of L,(x) drawn from the experimental 
values according to Fig. 2. This is fitted to the 
curves of Fig. 3 and the corresponding value of 
L,(0) obtained. y, can then be found from the 
relation 


L, (0) = 20 log coth y,. (13) 


In practice of course the value of L,(«) is small 
and for the usual case negligible and the measured 
curve L(x) is fitted directly to the curves in Fig.3 
and L,(0) and y, determined as before. 


1 It is, of course, permissible to determine L,(x) in Fig. 4 
by making a rough estimate for «x +, or taking it as zero, 
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4, The envelope function for highly absorbent 
samples 


(Case B) 


It is not practicable to use the method of the 
last section when L,(x) becomes comparable with 
L,(x); this occurs or can occur sometimes in the 
region when the value of y, exceeds 0.6 and cor- 
responds to values of L,(0) of about 5 dB or less. 
It is then better to use L(x) directly and to use 
x/k as the parameter. The envelope function L (x) 
is plotted against «x + y, for different values of 
o/k in Fig. 5.2 From this, it is possible to deter- 
mine (xx+y,) for the particular value «/k and 
the measured value of L(x) at distance x. y, can 
then be obtained by subtraction. 

It can be seen from Fig. 5 that there is theo- 
retically an upper limit of y, that can be meas- 
ured with a particular tube; the limit depends 
on the parameter «/k but is in general so near to 
100 per cent absorption coefficient as to be un- 
important. 


5. Measurements of the 
standing wave ratio 


The impedance tube used 
in the measurements describ- 
ed here is illustrated in Fig. 
6 and 7. The tube is made of 
brass, square in cross-section, 
and has an internal dimension 
of 7.5 cm, a wall thickness of 
0.9 em and alength of 170 cm. 
The ends of the tube are 
machined accurately square 
to the tube axis. A small 
probe tube 4.7 mm in diame- 
ter rests on the bottom of the 
main tube and is bent up- 
wards at the end so that the orifice is on the axis 
of the tube. The probe tube area is about 0.5 per 
cent of the main tube area and errors due to the 
probe can be expected to be smaller than this [13]. 
The impedance tube was terminated as required 
by a heavy steel plate (4 em thick), a brass plate 
(1.0 cm thick) or a sample contained in a holder; 
these plates and the sample holder can be seen 
in Fig. 7. Precautions were taken throughout the 
design to avoid unwanted vibrations in the micro- 
phone or the tube, to shield the microphone and 
to mount it resiliently. The loudspeaker box was 


2 Although for the measurement of the standing wave 
ratio of highly absorbent samples a recorder or meter with 
linear sensitivity is normally and preferably used instead 
of a high-speed level recorder with logarithmic sensitivity, 
the scale for L(x) in Fig. 5 has been given in dB for reason 
of uniformity. 
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Fig. 5. The envelope function L(x) plotted against «x + yp, 
for various values of «/k. 


>. 


Fig. 7. The end of the main tube with the different ter- 
minals and the sample holder. 
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decoupled from the main tube and from the probe 
tube where it passed through the speaker box. 
With these precautions and the use of a tuned 
amplifier (having a Q of about 6) the noise level 
when the orifice of the probe tube was sealed was 
almost 65 dB below the level at the maxima 
(see Fig. 10c). 

The attenuation constant « of the impedance 
tube and the parameter «/k are plotted against 
frequency in Fig. 8, and since «/k does not exceed 
0.005 at any point, the correction L,(x) is always 
negligible. The measured envelope curves from 
the examples of Fig. 9a and 9f can therefore be 
directly compared with the curves of Fig. 3; such 
comparison showed good agreement of experi- 
mental and theoretical values along the whole 
length of the tube. This provides an excellent 
check that the conditions in the tube are satis- 
factory. 


10:10°* 


—& Neper-cm7] 
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Fig. 8. Attenuation constant « of the main tube and ratio 
of x to the wave number k plotted against frequency. 


The usefulness of this method of checking is 
shown by results obtained when errors were 
deliberately introduced; these are illustrated in 
Fig. 10. The curve of Fig. 10a would normally 
be considered acceptable but the use of the en- 
velope curves (Fig. 3) shows that a mistake is 
being made. It was, in fact, caused by a deliberate 
increase in the speed of movement of the probe 
tube and a decrease in writing speed of the level 
recorder. It is interesting to observe that 44 dB 
would be the value of L,(0) used for calculating 
y, based on previous methods whereas 60 dB is 
the correct figure. During early work, the curve 
of Fig. 10b was obtained at low frequency and 
was recognized as incorrect since the curve passes 
the L,(x) line for complete reflection. Examina- 
tion of the output at the minimum showed a 
sinusoidal voltage with low noise. The curve was 
recorded again with a higher valued potentio- 
meter in the level recorder, the probe was sealed 
and the noise measured. The curves are given in 
Fig. 10c. The noise level was obviously low 
enough. Finally it was found that the trouble was 
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Fig. 9a--f. Records with envelope curves for a rigid ter- 
minal and different frequencies. 


caused by a very small coupling between the 
loudspeaker box and the main tube. After this 
had been removed, the correct result was ob- 
tained as shown in Fig. 10d. In most impedance 
tubes wherein high precision is sought, the check- 
ing necessary to ensure correct working takes a 
great deal of time and care. This can be greatly 
reduced by adopting the envelope technique de- 
scribed here. 

A further advantage of the use of the envelope 
curve over the whole tube is illustrated in Fig. 10e. 
Herein the minima of the recorded curve do not 
lie accurately on a smooth curve (undoubtedly 
owing to inherent small irregularities) but the use 
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of the family of envelope curves makes the task 
of choosing the envelope of best fit and the deter- 
mination of the best value of L,(0) much simpler 
and more accurate. 
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Fig. 10a---e. Records to illustrate the use of the envelope 
curves for detecting disturbances of the 
equipment. 


6. Measurement with a highly reflecting terminal 


In order to check the working conditions of an 
impedance tube or to determine the limits of 
accuracy for the difficult measurements of hard 
materials it is necessary or desirable to study the 
apparent impedance of a very hard terminal. The 
measured values of L,(0) for a ground steel plate 
4 cm thick as obtained with the equipment de- 
scribed above, are plotted against the frequency 
in Fig. 11, the technique of the envelope function 
being used in the measurement (cf. Fig. 9a---f). 
Several methods of fixing and sealing the plate 
were used and all gave essentially the same result. 
The corresponding absorption coefficients (cf. 
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Fig. 11. Values of L,(0), y, and absorption coefficient a, for 
a rigid terminal plotted against frequency. 
(a) average of measured values, 
(b) theoretical limit caused by heat conductivity 
only. 


curve a) are ().2 per cent at 100 c/s, 0.5 per cent 
at 500 c/s and 0.9 per cent at 2000 c/s. 

The theoretical limit is for this measurement 
plotted in Fig. 11 (curve b) also; it has been ob- 
tained by the method of Cremer [14] for the 
ideal case where the sound absorption is caused 
only by heat conductivity in a thin layer close to 
the surface of the metal. The slope of the two 
curves is the same but a difference of 15 dB in 
level indicates, as can be expected of course, that 
effects other than the heat conductivity deter- 
mine the measured limit. These may be the so- 
called edge effects at the terminal, disturbance 
of the sound field by the end of the probe tube or 
some unexplored effect of the equipment. 

In the first place the resolving power of the 
tube in connection with the possible error must 
be discussed. With a highly reflecting terminal, 
the first minimum occurs at a distance close to a 
quarter of a wavelength from the sample surface. 
For a given sound level at maxima, the level at 
the first pressure minimum is determined partly 
by the attenuation (x) of the tube and partly by 
sound absorption (y,) of the terminal. A factor of 
the resolving power D, may be defined as the 
ratio of the pressure at the first minimum with 
zero tube attenuation to the pressure at the first 
minimum for zero energy absorption at the ter- 
minal. 

Now for highly reflecting samples 


es oy | P,| Ric as —4y, ~<a 
wal (aa) ed, (14) 


and from eq. (10) 


(iis SMS fe (;) 


ah ah Ila, (15) 
and consequently for a constant a, the factor D, 
is proportional to the square root of the frequency. 


If, alternatively, L,(0) is determined from a 
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minimum other than the first, the corresponding 
resolving power D(x) becomes 


For constant a, and x, the factor D(x) is inversely 
proportional to the square root of the frequency. 
The use of minima at some distance from the 
sample is therefore a disadvantage at higher fre- 
quencies; however, the proposal made earlier 
that the record of the standing wave along the 
whole tube be used in determining L,(0) has the 
advantage that hidden disturbances can be re- 
cognized and eliminated. To avoid an over- 
estimate of the accuracy, consider the resolving 
power of the tube at a point about half its length. 

The measured absorption coefficient of Fig.11 
increases with frequency in such a way that the 
factor D(x) is independent of frequency. There- 
fore, the factor of the resolving power taken at 
the middle of the tube (x = 75 cm) is after eq. (16) 
and Fig. 8 and 11 in this case 


D(x = 75 cm) ~ 0.08. (17) 


If the factor D(x) is less than unity, it is obvious 
that the error of the measured absorption coef- 
ficient ay increases. Call the error of the absorp- 
tion coefficient ¢, and the error in the pressure 
ratio maximum to minimum ém. 


Then it yields 


(Pe oac” (Pes) Gee) 08) 


and Qo measured = A (1 ala Ea) (19) 
and éq and ¢é» are related 

= {1+ : 20 

£a = ( De) Em. ( ) 


For measurements of the envelope curve from 
a single tube record, the error may be estimated 
at +0.5 dB corresponding to ém=0.05. Repeti- 
tion and good experimental conditions may 
reduce these to + 0.1 dB, ¢n=0.01. The scatter- 
ing of measurements shown in Fig. 11 can be 
explained from eq.(17) and (20) and the tube 
conditions for determining the standing wave 
ratio with highly reflecting samples are as good as 
can be expected. 


7. The shift of position of the pressure minimum 


The complex value of the acoustic impedance 
is determined from the relative position of the 
pressure minimum and the sample as well as from 
the standing wave ratio (cf. eq. (4)). Normally 
the position d, of the first minimum (ef. Fig. 1) 


(16) 
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is chosen. The value of wy, is then found by the 
use of eq.(2); the wave length is preferably deter- 
mined from the average distance between two 


successive pressure minima (cs. b = seh: 

2 N-1 
The distance of a minimum depends to some 
extent on tube attenuation (eq.(2)) as discussed 
by Scott [5]. Let the change in position of the 
minimum with and without attenuation be 6x, 
then 


dx = — 575 sinh 2 (ax + y,) (21) 
and in all practical cases 6x is small for y, small. 
The effect of the attenuation on the position of 
minima is therefore small for hard materials but 
it is necessary to check and investigate the ac- 
curacy of the position of the first minimum for a 
highly reflecting sample. 

As the size and plane of the orifice of the probe 
tube vary frequently with different impedance 
tubes, four types of probe tube orifice have been 
used with the equipment described; their sections 
are given in Fig. 12. The difference 4x between 
the measured position d, and a quarter of a wave- 


length is plotted against the frequency in Fig. 13 
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Fig. 12. Cross-sections of probe tube orifices used in the 
measurements of Fig. 13. 
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Reproducibility of the order of +0.01 cm was 
achieved, and there was no significant difference 
by using different methods of fixing and sealing 
the end plate (steel 4 cm thick). 

The curves of Fig. 13 indicate that appreciable 
errors can be introduced by having the orifice 
perpendicular to the plane of the main tube (end 
pieces A and B) and that the effect is greater with 
larger orifices. End pieces C and D having orifices 


O00 200 
——=> Frequency [cls] 


400 600 1000 2000 


Fig. 13. Measured difference Ax of the position of the first 
minimum and a quarter of a wave length plotted 
against frequency for different tube orifices. 


parallel to the axis of the tube, gave values of Ax 
much smaller than A and B which moreover 
were nearly equal.* The end piece C was therefore 
chosen as the most suitable and was the one used 
elsewhere in the work. Comparison of Ax as meas- 
ured and 6x as calculated from eq. (21) show that 
6x can invariably be neglected. For example, at 
200 c/s, Ax is 0.09 cm and 6x is — 0.0005 em. 

An examination of the measurable values of 
w. can be made in a manner similar to the dis- 
cussion of the limits of y, given earlier. A practic- 
able experimental check of the condition of an 
impedance tube can be made by determining the 
frequency characteristics of L,(0) or a, and of 
Ax (as has been shown in Fig. 11 and 13). Such 
characteristics will provide a good specification 
of the accuracy possible from the particular im- 
pedance tube. 


8. Conclusions 


The envelope curve of the pressure-distance 
relation in an acoustic impedance tube as recorded 
with a high-speed level recorder can be deter- 
mined theoretically and its use is likely to im- 
prove the accuracy of the measurements. For 


3 It may be observed that the measured increase of Ax 
with lower frequencies does not necessarily mean that it is 
caused by the orifice of the probe tube and could not be 
avoided with other equipment. 
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highly reflecting samples, the family of curves for 
different envelopes is represented by taking the 
absorption of the sample (L,(0) or y,) as para- 
meter and for highly absorbent samples the 
curves have the ratio of attenuation to wave 
number (x/k) as parameter. Disturbing effects in 
the equipment can be rapidly detected by closing 
the tube with a rigid terminal, recording a 
pressure-distance curve and comparing this with 
theoretically derived curves. 

A definition of the resolving power of an im- 
pedance tube is proposed and discussed. The size 
and plane of the probe tube orifice is important — 
in particular the plane should be parallel to the 
tube axis. Test measurements with a rigid ter- 
minal whose results are plotted against frequency 
as shown in some examples provide a practicable 
method for specifying the tube performance and 
the accuracy of measurements possible with an 
impedance tube. 
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DIE MESSUNG DER TRITTSCHALLDAMMUNG VON DECKEN 


MIT SINUSFORMIGER ERREGUNG 


Von THOMAS LANGE 
IIL. Physikalisches Institut der Universitat Gottingen 


Zusammenfassung 

Es wird die Trittschalldammung von FuBbéden bzw. Decken in Gebauden gemessen. Die 
Anregung der Decke zu Schwingungen erfolgt sowohl mit dem iiblichen Hammerwerk als auch 
sinusformig mit einem elekiromagnetischen System. Die Differenz beider Trittschallpegel-Kurven 
stimmt mit dem Schlagspektrum des Hammerwerkes iiberein. Messungen auf dem Bau beweisen 
die Brauchbarkeit des Verfahrens zur Messung der Trittschalldammung von Decken bei sinus- 
férmiger Anregung, welches eine Reihe von Vorteilen bietet. Die Amplituden der Biegeschwin- 
gungen eimer gemessenen Decke nehmen bei punktweiser Erregung mit der Wurzel aus dem 
Abstand ab. 


Summary 


This paper describes investigations of the sound pressure level of impact sound. The vibrations 
of the ceilmg are produced both by the usual hammer machine and by an electromagnetic system 
for sinusoidal excitation. The difference of both sound pressure levels is in agreement with the 
spectrum of the impacts with the hammer. Measurements in buildings demonstrate the adequacy 
of the method of simusoidal excitation. This method has several advantages. In the case of excita- 
tion at one point the amplitudes of flexural waves on a floor decrease with the square root of 


the distance. 


Sommaire 


Mesure de la transmission sonore des bruits de choe dans les batiments. Les vibrations d’un 
plancher sont produites non seulement par la machine 4 frapper (& marteaux) mais aussi par un 
systéme électromagnétique excitant seulement une fréquence. La différence des niveaux sonores 
moyens produits dans le local de réception sexplique par le spectre des chocs des marteaux. 
Dans les batiments, des essais avec le systéme électromagnétique montrent lutilité de cette 
méthode d’excitation, qui posséde plusieurs avantages. Les amplitudes mesurées des vibrations 


de flexion d'un plancher diminuent avec la racine carrée de la distance du point excité. 


1. Einleitung _ 
Die Trittschalldammung von FuBbéden bzw. 


Decken in Gebauden wird nach einem genormten 
Verfahren [1]. [2] gemessen, indem man ein be- 
stimmtes Hammerwerk auf den FuBboden stellt. 
das stoBweise die Decke zu Schwingungen anregt. 
In dem darunter liegenden Raume mift man den 
mittleren absoluten Schallpegel in Frequenz- 
bandern von Oktavbreite. Im folgenden wird 
untersucht, wie die stoBweise Erregung durch 
eine sinusférmige ersetzt werden kann. 


2. Messungen mit dem Hammerwerk nach 
DIN 52210 


Zunachst sei tber einige Untersuchungen mit 
dem Hammerwerk berichtet, das den internatio- 
nalen Normen [2] entsprach bis auf folgende 
Punkte: Die Hammer bestanden aus Eisen. Jeder 
Hammer wurde etwa 250 ms nach dem Schlag 
wieder emporgehoben. Es interessiert dabei der 
Druck-Zeit-Verlauf wahrend eines Hammer- 


schlages: er wurde aufgenommen, indem der 
Hammer auf ein 2 em dickes, zylindrisches Eisen- 
stiick von 2 em Durchmesser schlug, welches den 
Druck auf einen darunter liegenden Piezokristall 
(X-Schnitt-Quarz mit Eigenfrequenzen ber 
0.3 MHz) iibertrug. Dieser war auf einem im Ver- 
gleich zum Hammer recht schweren Eisenklotz 
befestigt, welcher zur Schwingungsdampfung in 
Sand [3] gebettet war. Die durch ein einmaliges 
Herunterfallen des Hammers erzeugte piezo- 
elektrische Spannung wurde verstarkt (obere 
Grenzfrequenz 1 MHz) und auf einem Oszillo- 
graphen in ihrem zeitlichen Verlauf sichtbar ge- 
macht. Abb. | zeigt, daB nach jedem Schlag der 
Hammer wieder abspringt, erneut herunterfallt 
und so fort. Dasselbe Bild bietet ein auf den 
Hammer als Beschleunigungsempfanger gesetzter 
Piezokristall (Biegeschwinger unterhalb seiner 
Eigenfrequenz) (Abb. 2). Es entsteht eine Schlag- 
folge, bevor der Hammer von dem Antriebs- 
mechanismus wieder emporgehoben wird. Beim 
ersten Aufprall hat der Hammer 78,5°% seiner 
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Fallenergie verloren, der erste Nachschlag erfolgt 
nach 84 ms. 

Das in der Unterlage erzeugte Schlagspektrum 
wurde mit dem oben genannten Druckempfanger 
aufgenommen, der an einen Hochtonanalysator 
[4] mit logarithmischem Pegelschreiber ange- 
schlossen war. Die erhaltenen Spektren haben bei 


pe si 
50 ms 


Abb. 1. Der Druck-Zeit-Verlauf beim Schlag auf den Eisen- 
klotz. (Ordinatenmafstab nicht linear, um die 
kleinen Nachschlage sichtbar zu machen.) 


ewe lees 
50 ms 


Abb. 2 


. Der Beschleunigungs-Zeit-Verlauf beim Schlag auf 
den Eisenklotz. (OrdinatenmaBstab nicht linear, 
um die kleinen Nachschlige sichtbar zu machen.) 


der genormten Schlagfrequenz von 10 s~! eine 
Liniendichte von 2 Hz. Aus technischen Griinden 
tritt in dem benutzten MeBverfahren eine gerin- 


gere Liniendichte auf. Das Spektrum des Schlages 
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auf den Eisenklotz (Abb. 3) ist bis 20 kHz mit 
einer Abweichung von +5 dB konstant. 

In dem Spektrum des Schlages auf eine 
1x 0,6 m? groBe und 28 mm dicke Holzplatte 
mit 1 mm Linoleumbelag (Abb. 4a) sind die Dif- 
ferenzen zwischen Maxima und Minima gréfer 
(15 dB). AuBerdem tritt ab 5 kHz ein allgemeiner 
Abfall nach hohen Frequenzen auf, der etwa 
1 dB/kHz betragt. Er stimmt qualitativ mit der 
Theorie [5] ttberein;-quantitative Ubereinstim- 
mung -war nicht zu erwarten, da nicht alle Voraus- 
setzungen der Theorie erfillt waren. Die experi- 
mentelle Relaxationszeit betragt 2-10~4s. Nach- 
schlage traten an der Holzplatte nicht auf. 
Abb. 4b zeigt das Spektrum eines Schlages auf 
den Ejisenklotz, wenn ein 2 mm dickes Stick 
Schwammgummi zwischen Hammer und Unter- 
lage liegt. Der Abfall betragt hier 6 dB/kHz, die 
Relaxationszeit ist gréBer als 10~? s. 

Mit dem Hammerwerk, einem Siemensschen 
Schallpegelmesser und Oktavsieb wurde der Tritt- 
schallpegel einer Decke nach den Vorschriften 
DIN 52210 gemessen. Man kann wahrend der 
Messung mit dem Hammerwerk auf dem Fuf- 
boden herumfahren, anstatt es an mindestens 
sechs verschiedenen Stellen auf den FuBboden 
klopfen zu lassen. Dadurch erreicht man die- 
selben Ergebnisse, eine schnellere Messung, eine 
bessere Mittelung sowie einen Uberblick iiber die 
durch stehende Wellen hervorgerufene Schwan- 
kungen. Abb. 5 zeigt einen Ausschnitt aus einem 
so hergestellten MeBstreifen bei einer festen 
Oktavsiebeinstellung. Anstelle des Oktavsiebes 
kann man auch ein schmaleres Filter benutzen. 

Abb. 6 zeigt den Trittschallpegel in Abhangig- 
keit von der Frequenz fiir eine 32 cm dicke Decke. 
Diese Decke, eine alte Konstruktion (in dem 
Hause BunsenstraBe 7 des III. Physikalischen 
Institutes), ist 507 x 315 em? groB und besteht 
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Abb. 3. Das Spektrum des Schlages auf den Eisenklotz. 
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aus Zementbeton, der zwischen Doppel-T-Trager 
gestampft ist, aus Zementestrich und aus Lino- 
leum. Der darunter befindliche Raum hat ein 
Volumen von 40,5 m® und die ebenfalls auf Abb.6 
(unten) gezeigte Nachhallzeit. Die Reduktion auf 
die genormte Absorptionsflache in diesem Raum 
betriige 1,9 dB, ist aber im folgenden nicht beriick- 
sichtigt, da es sich nur um Vergleichsmessungen 
an ein und derselben Decke handelt. Die ge- 
strichelte Kurve auf Abb. 6 zeigt den Trittschall- 
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pegel derselben Decke, aber gemessen mit einem 
Filter von 50 Hz Halbwertsbreite (und 100 Hz 
Bandbreite bei einem Pegel 40 dB unter dem 
Maximum) [4]. Hier sind die MeBwerte bereits 
mit einer Korrektur versehen, welche die kon- 
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Abb. 4. Das Spektrum der Schlige auf die Holzplatte (a) 
und auf den Eisenklotz mit Gummizwischen- 


schicht (b). 
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Abb. 5. Schwankungen des Trittschallpegels, wenn man die 
Aufstellung des Hammerwerkes verschiebt. 
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Abb. 6. Der Trittschallpegel, ausgezogene Kurve bei Oktav- 
sieb-Analyse; gestrichelt bei Analyse mit 50-Hz- 
Filter, umgerechnet auf Oktavbandbreite. 


stante Bandbreite umrechnet in eine mit der Fre- 
quenz wachsende, wie sie ein Oktavsieb besitzt. 
Dann stimmen beide Kurven innerhalb der MeB- 
genauigkeit iiberein. Soweit die Messungen nach 
der Hammerwerk-Methode. 


3. Verfahren zur Messung mit sinusformiger 
Anregung 

Zur sinusférmigen Anregung von Kérperschall 
wurde ein elektromagnetisches System gebaut. 
Es besteht aus einem 13 kg schweren Permanent- 
magneten (Type Nr. 20-8617a der Deutschen 
Edelstahl-Werke, Dortmund), der in einem ring- 
formigen Spalt von 62 mm Durchmesser, 1,6 mm 
Breite und 10 mm Tiefe eine KraftfluBdichte von 
etwa 1,3 V-s-m ® besitzt. In diesem Spalt sitzt 
eine Spule, die auf einen zylinderférmigen Spulen- 
trager aus Messing gewickelt ist und durch zwei 
Blattfedern so gehalten wird, da sie in dem 
Magnetspalt frei schwingen kann. Abb. 7 zeigt 
das System, das also einem permanent-dynami- 
schen Lautsprecher ahnlich aufgebaut ist. Nur ist 
anstelle der Membran ein starrer Fu angebracht 
(auswechselbar mit Durchmessern von 30 bis 
130 mm), auf den der ganze Magnet gestellt we- 
den kann. Er wird dann von den Blattfedern ge- 
tragen. Diese wurden so berechnet, daB sie zu- 
sammen mit der Masse des Magneten eine Eigen- 
frequenz von 15,9 Hz ergeben. Die tatsachliche 
Eigenfrequenz liegt bei 15,6 Hz. Diese Uberein- 
stimmung ist iiber Erwarten gut. 
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Abb. 7. Der elektromagnetische Kérperschallsender schrig 
von unten (etwa 1/, natiirlicher GréBe). 


Wird dieser Kérperschallsender auf eine Unter- 
lage, z. B. den FuBboden, gestellt und durch seine 
Spule ein Wechselstrom I aus einem Verstarker 
geschickt, so wirkt zwischen Sender und Unter- 
lage eine Kraft K,die dem Strom proportional ist: 
K=alI. Da « hier nicht von der Frequenz ab- 
hangt, ist der Effektivwert der Kraft leicht fre- 
quenzunabhangig zu machen. (Dies gilt auch bei 
und unterhalb der Eigenfrequenz. Da in der 
Eigenfrequenz jedoch besonders grobe Ampli- 
tuden auftreten, werden die Verluste gréfer. 
AuBerdem ist die Justierung der Schwingspule 
fiir groBe Amplituden sehr schwierig. Es emp- 
fiehlt sich daher, den Sender nur oberhalb seiner 
Eigenfrequenz zu benutzen.) 

Die Frequenzunabhangigkeit des Effektiv- 
wertes der Kraft zwischen Sender und Unterlage 
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Abb. 8. Die Frequenzkurve des gesamten Systems. 
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wurde gemessen, indem der Sender auf einen 
piezoelektrischen Druckempfanger gestellt wurde. 
Das Ergebnis zeigt die Abb. 8. Die untere Grenz- 
frequenz von 75 Hz ist durch den Kraftverstarker 
bedingt, die obere bei 1,5 kHz durch die Reso- 
nanz des piezoelektrischen Empfangers, der aus 
Seignettesalz-Kristallen zwischen zwei als Masse 
wirkenden Messingplatten besteht; durch das 
Aufsetzen des Kérperschallsenders wurde die 
Masse vergroRert und seine Eigenfrequenz von 
3,5 auf 1,8 kHz gesenkt. Ein stérender EinfluB 
von Blattfeder-Resonanzen wurde nicht beobach- 
tet. Durch Verwendung eines Verstarkers mit 
besserem Ausgangsiibertrager konnte die Fre- 
quenzkurve noch verbessert werden. 

Die Proportionalitatskonstante « zwischen 
Kraft und Strom des Kérperschallsenders wurde 
wie tiblich in einem statischen Versuch bestimmt. 
Man 1aBt zwischen Spule und Magnet eine be- 
kannte Gleichkraft wirken, indem man ein Ge- 
wicht auf den Magneten setzt. Die hierdurch ver- 
ursachte Stérung der Gleichgewichtslage zwischen 
Spule und Magnet wird durch einen Gleichstrom 
kompensiert. Der so erhaltene Wert fiir « lautet 


x= 14,1 Newton/A=1,44 kp/A +4%. 


Urspriinglich war « gréBer und sank im Laufe 
der Zeit infolge Alterung des Magneten auf den 
oben angegebenen Wert. 

Oberhalb 1 kHz strahlte der Magnet merklich 
auch Luftschall ab, den man durch eine tiber den 
Magneten gestiilpte Dammhaube geniigend ver- 
mindern kann. Die Dammbhaube besteht aus 
einem doppelwandigen GefaB, das mit Sand in 
einer 22 mm dicken Schicht gefillt ist. 

Die zur Messung des Trittschallpegels benutzate 
Apparatur besteht auf der Sendeseite aus einem 
Schwebungssummer mit Verstarker (90 W) und 
dem beschriebenen elektromagnetischen Anrege- 
system. Die Empfangsapparatur setzt sich aus 
einem DIN-Schallpegelmesser (Siemens & Halske 
AG) und einem daran angeschlossenen (Neu- 
mannschen) Pegelschreiber zusammen. Die An- 
triebsachse des Pegelschreibers wird mit der 
Achse des Drehkondensators im Schwebungs- 
summer starr verbunden. Ohne Frequenzmodula- 
tion (Wobbelzusatz) des Summers zeichnet der 
Pegelschreiber die stehenden Wellen des Luft- 
schalles im Empfangsraum iiber der Frequenz 
auf, ahnlich wie bei der Messung der Luftschall- 
dammung. Nun wird der Magnet auf eine andere 
Stelle auf den FuBboden gesetzt, das Registrier- 
papier zuriickgerollt und eine weitere Kurve iiber 
die erste geschrieben. Abb. 9 zeigt, daB sich beide 
Kurven bis auf die stehenden Wellen decken. Auf 
diese Weise werden sechs Kurven iibereinander 


oe 


ACUSTICA 
Vol. 3 (1953) 


a ete 


oo ea 


TH. LANGE: MESSUNG DER TRITTSCHALLDAMMUNG VON DECKEN 


165 


Bog 


Abb. 9. Ein Registrierstreifen des Trittschallpegels bei sinusférmiger Anregung. Zwei Kurven, gewonnen 


an verschiedenen Anregungsstellen. 


gezeichnet, deren Mittelwert in Abhangigkeit von 
der Frequenz graphisch ermittelt wird. 

An der oben beschriebenen Decke wurde der 
Trittschallpegel gemessen; die Anregung erfolgte 
durch das elektromagnetische System mit ver- 
schiedenen FuBdurchmessern von 30 bis 130 mm. 
Die Dicke der als Fu dienenden Dural-Scheiben 
betrug 15 mm. Eine Abhangigkeit von der An- 
regungsflache ist zu erwarten, sobald deren 
Durchmesser mit der Biegewellenlange der Decke 
vergleichbar wird, sofern der Fuf als starr gegen 
Biegeschwingungen zu betrachten ist. Dies ist 
aber nicht erwiinscht und eine derartige Starrheit 
1aBt sich wohl gar nicht mit Konstruktionen er- 
reichen, deren Flachengewicht soviel kleiner ist 


als das der Decke. Die MeBergebnisse (Abb. 10) 
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Abb. 10. Der Trittschallpegel bei sinusférmiger Anregung. 
Sechs verschieden groBe kreisférmige Anregungs- 
fliachen mit Durchmessern von 3 bis 13 cm. 


zeigen keine Abhangigkeit des gemessenen Tritt- 
schallpegels vom Fubdurchmesser. Vielmehr wird 
die Streugrenze und damit die MeBgenauigkeit 
aus dieser Darstellung deutlich. Sie betragt 
+2 dB. 


4, Vergleich der Trittschallmessungen bei schlag- 
artiger und bei sinusformiger Anregung 


Ersetzt man die Anregung einer Decke zu 
Biegeschwingungen mit dem Hammerwerk durch 
eine sinusférmige Anregung mit bekannter Kraft 
K, soist der Trittschallpegel L auf die GréBe dieser 
Kraft zu beziehen. Man erhalt so mit der tiblichen 
Reduktion auf konstante Absorptionsflache A, im 
Empfangsraum im logarithmischen Ma’: 


Lin = 20 log p/p y —20 log K/K, + 10 log A/Ap. 

Als Bezugswerte wurden die GréBen 

Po = 2: 10—> Newton/m? = 2-107* ubar = 

rl OG dymicm* 

K, = 10 Newton = 10° dyn = 1,02 kp 

Fl, == WO yane 
gewahlt, weil sie Werte in derselben GréSenord- 
nung ergeben wie bei Anregung mit dem genorm- 
ten Hammerwerk. Man kénnte auch andere Be- 
zugswerte festlegen. 

Die Definition des Norm-Trittschallpegels, das 
heiBt bei der Messung mit dem Hammerwerk, 
lautet (in der Schreibweise von KosTEN u. a. [9]): 
LNorm (M = A,/A) = 20 log p/p, + 10 log A/Ay 
Po = 2°10-> Newton/m? = 2-104 pbar = 

= 2. 10—* dynj/em?, 
Ag = 10im?. 
Die Reduktion auf eine bekannte und frequenz- 
konstante Kraft entfallt hier, da man sich auf die 
durch das Hammerwerk gegebene Kraft mit 
seiner durch die értliche Elastizitat bedingten 
Frequenzabhangigkeit bezieht. 


An der oben schon beschriebenen Decke wurde 
ein Vergleich der Trittschallpegel bei Anregung 


mit 
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Abb. 11. Der Trittschallpegel; ausgezogene Kurve bei An- 
regung mit Hammerwerk; gestrichelt bei sinus- 
formiger Anregung; Kreuze bei sinusférmiger An- 
regung, aber umgerechnet auf Oktavbandbreite. 


mit dem Hammerwerk und mit dem Magnet- 
system durchgefiihrt. Abb. 11 zeigt glatt ausge- 
zogen den Trittschallpegel bei stoBweiser Anre- 
gung durch das Hammerwerk und gestrichelt bei 
sinusférmiger Anregung durch das Magnetsystem. 
Beide Kurven sind zunachst noch nicht mitein- 
ander vergleichbar, da bei der Oktavsiebanalyse 
die Filterbreite proportional der Frequenz wachst, 
wahrend sie bei der Sinuserregung konstant 
bleibt. Die mit dem Magnetsystem gewonnene 
Kurve wurde daher von konstanter Filterbreite 
auf Oktavbreite umgerechnet. Das Ergebnis ist 
ebenfalls in Abb. 11 als ausgezogene Kurve mit 
Kreuzen eingetragen. Sie stimmt noch nicht mit 
dem Trittschallpegel bei Anregung mit dem 
Hammerwerk iiberein, da die Anregungsspektren 
beide Male verschieden sind. 

Wenn dies so ist, sollte die Differenz der Kur- 
ven des Trittschallpegels, einmal gemessen mit 
dem Hammerwerk, das andere Mal mit sinus- 
formiger Anregung mit dem Magneten und auf 
Oktavbandbreite bezogen, das Anregungsspek- 
trum des Hammerwerkes darstellen. Diese Diffe- 
renz ist als ausgezogene Kurve in Abb. 12 darge- 
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Abb. 12. Das Spektrum des Hammerschlages auf die Decke. 
Ausgezogene Kurve: Differenz der Trittschall- 
kurven bei Anregung mit Magnetsystem bzw. 
Hammerwerk; Punkte aus Messung mit unter- 
gelegtem Piezoquarz. 
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stellt und wird zur Bestatigung, da keine weite- 
ren Ejinfliisse von wesentlicher Bedeutung da sind, 
mit dem direkt gemessenen Spektrum verglichen. 
Letzteres wurde, wie bereits beschrieben, mit 
einem unter den Hammer auf die Decke gelegten 
Piezoquarz gemessen. Die Mefwerte sind als 
Punkte in Abb. 12 eingezeichnet. Man sieht, daB 
der Frequenzverlauf des Hammerspektrums pro 
Oktave derselbe ist, wie er aus der Differenz der 
Trittschallpegel-Kurven. gewonnen wurde. Die 
Streuung der MeSpunkte riihrt davon her, da 
jeder eingezeichnete MeSpunkt nur einem ein- 
zigen Hammerschlag entspricht, wahrend die 
ausgezogene Kurve aus dem Mittelwert sehr 
vieler Schlage gewonnen wurde. 


5. Messungen auf dem Bau 


Die Brauchbarkeit des geschilderten Verfahrens 
zur Messung des Trittschallpegels bei sinusférmi- 
ger Anregung der Decke mit dem Magnetsystem 
wurde durch Messungen auf dem Bau bestatigt. 
Zur Verfiigung stand eine sogenannte Phénix- 
deckenkonstruktion, die im wesentlichen aus 
Stahlbetonrippen besteht, deren Zwischenraume 
zwischen den Rippen mit Hiittenbims ausgefiillt 
sind. Das Flachengewicht dieser Konstruktion 
betragt 426 kg/m?. Es wurden Messungen an 
zwei gleichen, nur in der GréBe verschiedenen 
Decken ausgefiihrt : S, = 19,6 m? und S, = 12,2 m?. 
Unter diesen Decken befanden sich Raume von 
V,=50,7 m3 und V,=31,5 m*? Volumen, deren 
Nachhallzeiten auf den Abb. 13 und 14 unten 
eingezeichnet sind. Diese Bilder zeigen auBerdem 
den gemessenen Trittschallpegel (ohne Reduktion 
auf eine bestimmte Nachhallzeit oder Absorp- 
tionsflache) als ausgezogene Kurven bei Anre- 
gung mit dem Hammerwerk und gestrichelt bei 
Anregung mit dem Magnetsystem, bezogen auf 
eine anregende Kraft von 10 Newton= 10° dyn. 
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Abb. 13. Messung im Bau: GroBes Zimmer. Ausgezogene 
Kurven: unten Nachhallzeit, oben Trittschall- 
kurve bei Anregung mit Hammerwerk; gestrichelt 
bei sinusférmiger Anregung. 
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Abb. 14. Messungen im Bau: Kleines Zimmer. Ausgezogene 
Kurven: unten Nachhallzeit, oben Trittschall- 
kurve bei Anregung mit Hammerwerk; gestrichelt 
bei sinusférmiger Anregung. 


Die Messung erfolgte mit einem Anregungsstrom 
von 2 A, durch die Schwingspule des Magneten, 
entsprechend einer Kraft von 24,2 Newton, so daB 
der empfangene Pegel durchweg tiber demjenigen 
lag, der mit dem Hammerwerk gewonnen wurde. 
Dies ist von Bedeutung fiir den Stérabstand 
(d.i. die Differenz von Nutz- und Stérpegel), der 
méglichst groB sein soll. 

Uber die Abhangigkeit der Trittschalldam- 
mung von der GriBe der Priifdecke ist bisher 
wenig bekannt. In theoretischen Betrachtungen 
wurden teils [6] der mathematischen Einfachheit 
halber ebene Wellen, teils [5], [7] Zylinderwellen 
angenommen. Experimentell wurde die auf einer 
Decke herrschende Wellenform ermittelt, indem 
die Decke mit dem Magnetsystem zu Schwin- 
gungen von 3+ 0,15 kHz angeregt wurde. Die 
Amplitude wurde in Abhangigkeit von der Ent- 
fernung vom Sender mittels eines piezoelektri- 
schen Beschleunigungsmessers [8] abgetastet. Der 
Sender stand in der Mitte des Zimmers. (Bei An- 
regung in der Nahe des Randes wird die Streuung 
der MeBpunkte noch gréBer.) Um beurteilen zu 
kénnen, mit welcher Potenz vom Abstand r die 
Amplitude abnimmt, ist die Amplitudenvertei- 
lung in Abb. 15 doppelt logarithmisch dargestellt. 
Die eingezeichnete Gerade entspricht einer Ab- 
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Abb. 15. Die Amplitudenverteilung auf der Deckendiago- 


nalen bei Anregung in der Mitte (0 cm) mit 
3+0,15 kHz. 
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nahme mit r—”, wie sie fiir Zylinderwellen gilt, 
und paBt recht gut zu den MeBpunkten. (Es han- 
delt sich um Messungen in einem Neubau auf der 
oben beschriebenen Phénix-Decke.) 


6. Diskussion 


Bei der Messung der Trittschalldammung mit 
sinusférmiger Anregung handelt es sich um eine 
physikalisch durchsichtige Methode zur Prifung 
der Eigenschaften von Decken. Sie unterscheidet 
sich von der ttblichen Messung mit einem ge- 
normten Hammerwerk wesentlich durch die Be- 
nutzung eins konstanten (oder doch wenigstens 
bekannten) Anregungsspektrums. 

Weitere Vorteile der stationdren Anregung 
sind: 1. ein groBer Stérpegelabstand; 2. die Még- 
lichkeit, die gleichen Apparate benutzen zu kén- 
nen, welche zur Messung der Luftschalld4mmung 
gebraucht werden (nur tritt an die Stelle des 
Lautsprechers ein Kérperschallsender, der iibri- 
gens einem Lautsprecher Ahnlich konstruiert 
werden kann); 3. die Méglichkeit der automati- 
schen Registrierung der MeBergebnisse, wie sie 
bei der Messung der Luftschalldammung vielfach 
ublich ist; 4. die Definition einer echten Isolation 
fiir den Trittschall. Diese Dammung miiBbte die 
Form besitzen (wieder in der Schreibweise von 
KOSTEN u.a. [9]): 


D(M =S/A) = 201lo us +10 log S/A. 
5 Sop i) 


Hierin bedeuten S die Deckenflache, K die auf 
die Decke an einer Stelle ausgeiibte Wechselkraft 
und p den Schalldruck im Empfangsraum mit der 
gesamten Absorptionsflache A. Entweder setzt 
man fiir S, die Flache S der untersuchten Decke 
ein. Dann betrachtet man den Mittelwert des 
anregenden Druckes iiber die ganze Deckenflache, 
wahrend dieser Druck in Wirklichkeit nur auf 
einer sehr kleinen Anregungsflache wirksam ist. 
Oder man setzt fiir S, einen konstanten Bezugs- 
wert ein, wie es oben in den Abschnitten 4 und 5 
geschehen ist. In beiden Fallen erhalt man als 
Dammzahl Werte, welche imstande sind, die 
Trittschalldammung einer Baukonstruktion zu 
charakterisieren. Diese Werte sind abhangig von 
der Frequenz, von den elastischen Kigenschaften 
der Decke (E-Modul, Biegesteifigkeit, innere 
Schwingungsdampfung), ihrer Masse, Dicke, Kin- 
spannung am Rande und von ihrer GréBe. 

Ich danke Herrn Prof. Dr. Erwin MEYER fiir 
die Anregung zu dieser Arbeit, sein férderndes 
Interesse daran und viele wertvolle Ratschlage. 

Fiir die Durchfithrung der Aufgabe stellte das 
Bundesministerium fiir Wohnungsbau die Mittel 


zur Verfiigung. (Eingegangen am 24. April 1952.) 
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THE ACCEPTABILITY OF SPEECH AND MUSIC 
WITH A SINGLE ARTIFICIAL ECHO 
by R.W. MUNCEY, A. F. B. NICKSON and P. DUBOUT 


Division of Building Research, Commonwealth Scientific and Industrial Research Organization, Australia 


Summary 

Groups of subjects seated in a room of 0.15 seconds reverberation time assessed whether or 
not a single artificial echo, added to speech or music previously heard with little or no echo, was 
disturbing. The artificial echo varied in intensity level from equal to the direct sound to —40 dB 
and in delay from 30 to 1000 milliseconds. The results obtained for speech differ from those of 
HAAS owing to the different reverberation time. An acceptable-echo-level/delay criterion deduced 
from the results is shown to be in agreement with pulse work in motion picture theatres. 


Sommaire 

On a demandé 4 des groupes de sujets assis dans une salle ayant un temps de réverbération 
de 0,15 seconde, d’estimer si un écho artificiel simple ajouté 4 de la parole ou de la musique 
écoutée auparavant sans écho ou avec un faible écho, génait ou non l’audition. On a fait varier 
Vintensité de l’écho artificiel entre ’intensité du son direct et 40 décibels au-dessous de cette 
intensité, et le retard de l’écho de 30 4 1000 millisecondes. Les résultats obtenus par la parole 
sont différents des résultats de Haas, car le temps de réverbération n’est pas le méme dans les 
deux cas. On déduit des résultats une valeur acceptable du rapport niveau/retard de l’écho, et 
on montre qu’elle est en bon accord avec les résultats des essais effectués avec des impulsions 
dans des cinémas. 


Zusammenfassung 


Gruppen von Versuchspersonen, die in einem Raum mit einer Nachhallzeit von 0,15 Sekunden 
saBen, ermittelten, ob ein Einfachecho, das zu Sprache oder Musik hinzugefiigt wurde, die vorher 
mit wenig oder keinem Echo gehért wurde, stérend war oder nicht. 

Das kiinstliche Echo wurde in der Schallstiirke im Vergleich zum direkten Schall von 0---—40dB 
und in der Verzégerungszeit von 30---1000 Millisekunden variiert. 

Die Ergebnisse fiir Sprache unterscheiden sich von denen von Haas wegen der unterschied- 
lichen Nachhallzeit. 

Ein angemessenes Kriterium fiir das Verhiltnis Echointensitaét zu Verzégerungszeit, das aus 
den Ergebnissen abgeleitet wurde, ist, wie gezeigt wird, in Ubereinstimmung mit Impuls- 
messungen in Filmtheatern. 


1. Introduction 


Sounds as heard inside a building consist of an 
initial direct sound followed by a succession of 
echoes of varying intensity and delay. Haas [1], 
[2] measured the intelligibility of speech when 
the initial sound was followed by one artificial 


echo and found that for each echo intensity level 
there occurred a “critical delay difference” at 
which 50 per cent of his listeners could hear a dis- 
turbing echo. The major part of his listening tests 
were done in a room of 0.8 s reverberation time. 
An increase in the delay caused an increase in the 
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number who found the echo disturbing. He found 
that the critical delay difference was inversely 
proportional to the number of syllables per second, 
was increased by decreasing the echo intensity 
level, was independent of the direction from 
which the echo came and was increased by an 
increase in reverberation time. 

These results have been used and extrapolated 
by Boutr and Doak [3] to provide a tentative 
criterion for the short-term (0 — 400 ms) transient 
response of an auditorium for speech. Since there 
is a good deal of interest also in concert halls it 
was decided to continue the pioneering work of 
Haas into the field of musical sounds and to 
check how this would affect the Bott and Doak 
criterion. 


2. Equipment 


The equipment has been designed to enable 
artificial echoes to be generated by the method 
used by Haas, i.e. by recording on a magnetic 
tape and varying the distance between record 
and replay heads to obtain the required delays. 
Since multiple artificial echoes are to be investi- 
gated later it was desirable to arrange a succession 
of replay heads. To obtain small variations in the 
echo time a tape speed of 1.5 m/s has been 
adopted and the distance between replay heads is 
continuously variable to give delays from 22 to 
1000 ms. The physical size of the heads is such 
that delays must be spaced at least 22 ms apart. 
The erase and the record amplifiers of an E.M.I. 
type B.T.R.I. tape recorder are used to record 
the signal on the tape, E.M.I. heads are used for 
replay and laboratory-made preamplifiers feed 
the control panel and thence Williamson am- 
plifiers and Australian Rola type 12-OX loud- 
speakers mounted in vented enclosures. The over- 
all performance of the equipment from 60 to 
8000 c/s can be considered adequate. The general 
layout is shown in Fig. 1. 

It was decided early in the work that a definite 
portion of music or speech should be repeated so 
that observers became familiar with it and could 
then make more rapid and consistent judgments. 
It is possible that by this a systematic error is 
introduced since a particular phrase may contain 
a critical combination, but this does not seem 
likely because the results obtained follow a de- 
finite and reasonable pattern. In the tests to be 
discussed an excerpt of fourteen seconds dura- 
tion was repeated for each type of sound, of 
which there were six: slow and fast speech, slow 
and fast string music and slow and fast organ 
music. 
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Each was recorded with a minimum of echoes. 
Speech was recorded standing close to a direc- 
tional microphone in the room used for the listen- 
ing tests, string music was recorded from a quar- 
tette ensemble playing out of doors and organ 
music was obtained with a directional micro- 
phone placed close to one speaker of a Compton 


| Listening room 


Operator's table 


Fig. 1. General layout. 


electronic organ. The tempo for the music can 
best be judged from the melodies given in Fig.2, 
each melody taking fourteen seconds. The ex- 
cerpts of speech were:—slow (about 2% syllables 
per second): ... “two light chairs had been over- 
turned, books and papers from the desk strewed 
the floor. The grandfather clock which should 
have stood sentinel on the left as one entered, 
had fallen”... and fast (about 5 syllables per 
second): ... “but for Henry Wilson, even though 


Fig. 2. Melodies of music used: 
(a) string music—fast tempo, 
(b) string music—slow tempo, 
(c) organ music—fast tempo, 
(d) organ music—slow tempo. 
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he had only just brought to the conclusion his 
most spectacular case, it was hardly likely that 
a single village eye would be batted for his arrival. 
Unless indeed a crime were to occur, but Michael 
had yet to see a place which less suggested 
crime’... 


3. Technique 


The subjects at each session were seated close 
together in a room 6x 4x 3 m® and the sound 
and the echo were supplied from two loud- 
speakers placed side by side. It is considered that 
Haas has satisfactorily demonstrated that the 
result with two such sources will be comparable 
with that for the arrival of the echo from any other 
direction. (For this test Haas was working out of 
doors.) The actual intensity level of the sound and 
its echo will of course vary slightly from seat to 
seat but this variation was shown to be less than 
2 dB and variability thereby introduced will be far 
less than the variation of judgment of the obser- 
vers. The subjects were asked in each test to decide 
whether they were disturbed by the echo (per- 
ceptible non-disturbing echoes were to be rated 
not-disturbed) and to record their judgment by 
moving a switch. Electrical means were used to 
obtain the total number recording ‘disturbed’ and 
this was displayed as a meter reading at the 
operator’s bench and recorded by him. Two sepa- 
rate teams were used as observers. The first team 
comprised the eight persons of the acoustics (and 
allied) group of the Division. Thirty persons taken 
from the remainder of the staff of the Division 
formed the second team from which twenty sub- 
jects were selected at random before any session. 
It was considered that a comparison of the results 
of the two groups would check that the learning 
period was about two sessions only and that more 
experienced subjects would not give a different 
result in the manner reported by HopxKinson [4] 
for glare experiments. 

At any one session, about eighteen repetitions 
of a particular sound type were presented with 
various combinations of echo delay and: level in 
a randomised order. This was repeated for three 
other sound types, again in a randomised order. 
The change of echo.was made with the sound 
muted and the new combination presented by 
opening the’shorting switch. The total sound level 
was arranged to be of constant intensity, i.e. it 
was assumed after HAAs that the ear added the 
sound levels incoherently. 

The sounds were presented at a level of 80 dB 
(peak) above 0.0002 dyne/cm? being that desired 
by the subjects for comfortable listening. Noise 
level in the room was about 35 dB. The subjects 
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Fig. 3. Reverberation time/frequency characteristics of lis- 
tening rooms; 
reverberant room, 
——-—-—-—— non-reverberant room. 
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500 ms 
Fig. 4. Response of listening rooms to 20 ms pulses of 

sound, 

(a) Non-reverberant room, microphone near loud- 
speaker; 

(b) non-reverberant room, microphone near centre 
of room; 

(c) non-reverberant room, microphone near rear 
wall; 

(d) reverberant room, microphone near centre of 
room. 
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were asked to record their votes during the period 
of silence between repetitions; voting during pre- 
sentation of a combination was discouraged. Ses-. 
sions were limited to 30 min to avoid fatigue 
effects. 

Finally, as a check on the difference between 
the results of this experiment and those of Haas, 
some experiments were made with a similar tech- 
nique in a room about 7 x 5 x 3 m? and of about 
1.0 s reverberation time. 
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two groups gave the same answer except that in 
the group of 20 a percentage of listeners was dis- 
turbed by sound with no echoes whereas the 
group of 8 voted not disturbed. For the fast sound 
types the same information is plotted in Fig. 7a---c 
in the form used by Bout and Doak giving vary- 
ing relations between echo intensity level and 
delays for 20, 50 and 80 per cent of subjects hear- 
ing a disturbing echo. It is considered that results 
for 20 and 80 per cent can be more accurately 
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Fig. 5. Percentages of listeners disturbed by single echoes of various sound types (non-reverberant room); 


(a) speech—fast tempo, 


(b) string music—fast tempo, 
(c) organ music—fast tempo, 


4. Results 


Reverberation-time/frequency curves for the 
two rooms used in the experiments are given in 
Fig. 3; pulse pictures obtained in the two rooms 
are reproduced in Fig. 4. 

The percentage of subjects registering a dis- 
turbing echo is plotted as a function of echo delay 
and echo intensity level for the various types of 
sound in Fig. 5a---f for the non-reverberant room 
and in Fig. 6a---e for the reverberant room. The 


(d) speech—slow tempo, 
(e) string music—slow tempo, 
(f) organ music—slow tempo. 


assessed than those for 10 and 90, as used by 
Bott and Doak. 

Each echo combination was presented at least 
four times except when a reading of all subjects 
disturbed was obtained at the first two or three 
presentations (further presentation was then con- 
sidered unnecessary). There was still a good deal 
of variation but with four repetitions the 95 per 
cent fiducial limits were about +10 per cent. 
This is adequate since factors like type of listener 
and tempo of the sound will undoubtedly intro- 
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Fig. 6. Percentages of listeners disturbed by single echoes 
of various sound types (reverberant room); 
(a) speech—fast tempo, 
(b) string music—fast tempo, 
(c) organ music—fast tempo. 
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the speech or music is slow and, as might be ex- 
pected from reverberation-time criteria, increas- 
ing echoes are tolerated as one moves along the 
range from speech to string and to organ music. 
Since obviously one does not design a hall for 
slow music or speech only, work was concentrated 
on fast speech and music and most of the dis- 
cussion is limited to that. 

Throughout all tests the sound was occasionally 
presented without echo.For speech in the non- 
reverberant room, subjects generally recorded 
not-disturbed, the average being about two per 
cent disturbed. For speech in the reverberant 
room and for music in both rooms a significant 
number still recorded a disturbing echo. The 
reverberation undoubtedly contributed but the 
psychological factor mentioned by Has is likely 
to be a prime factor, particularly for music: “the 
observer thought that almost always an echo was 
transmitted and he endeavoured therefore to 
avoid the blame of being inattentive” —(EHREN- 
BERG’s translation). Possibly also for music the 
ear is used to echoes and therefore prefers them, 
so that lack of echoes gives an unusual effect rated 
by some as disturbing. Subjects complained of 
the difficulty of detecting echoes, particularly 
with the slow organ music and the percentage 
disturbed when no echo is added is a rough meas- 
ure of the difficulty of detecting echoes. It seems 
reasonable then to treat the percentage disturbed 
with no echo, obtained from the group of 20, as 
the line of zero disturbance (since in a hall people 
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Fig. 7. Acceptable echoes for various percentages of listeners disturbed (non-reverberant room); 
(a) speech—fast tempo, (b) string music—fast tempo, (c) organ music—fast tempo. 


duce errors equally large. The various average 
curves as drawn are a reasonable mean and fall 
inside the fiducial limits. 


5. Discussion 


The acceptability of an echo added to a sound 
originally without echo varies with its intensity 
and delay and with the type of sound and its 
tempo. An increase in intensity level or delay 
tends to decrease the acceptability. Louder echoes 
can be added without listener-disturbance when 


are less conscious of and are not actually expect- 
ing echoes) and to rate percentages on this basis. 
This has been done on the right-hand scale in 
Fig. 5 and 6 and these percentages are used in 
Fig. 7. 

The results do not agree exactly with Haas’ 
claims of the acceptability of short-delay echoes 
of the order of + 10 dB although + 10 dB echoes 
are almost as acceptable as those of the + 5, 0 and 
—5 dB values. This is probably partly due to 
the psychological factor mentioned earlier coupled 


— 
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with subtle changes in the reproduction quality. 
However from the work of PARKIN and his as- 
sociates [5], [6] it is quite obvious that sound re- 
inforcement of the order of +10 dB will pass 
unnoticed when the audience is unconscious of 
its use. When a speech reinforcement system is 
in use the optical sense is of importance since this 
sense will demand that the source is other than 
the reinforcement speaker and may well explain 
a 10 dB difference in effect. 

Haas’ results concerning the echo disturbance 
as a function of a given delay and intensity show 
less disturbance than do those obtained in the 
non-reverberant room. The lower reverberation 
is the obvious reason for this fact. Haas noticed 
that a decrease in reverberation time decreased 
his “critical delay time” for echoes of intensity 
level equal to the original sound. The effect is far 
more spectacular however for echoes of the order 
of —10 dB. Check measurements made in a room 
of reverberation time about one second gave 
results comparable with those of Haas. It seems 
likely that in such a room the echoes from the 
bounding surfaces are louder than the added echo 
and being more natural are acceptable. 

The new measurements require for speech a 
more stringent criterion than that previously pro- 
posed. P.H. Parkin of the Building Research 
Station, Great Britain has informed the authors 
that the Bor and Doak criterion agrees with 
subjective results at the Shakespeare Memorial 
Theatre, Stratford-on-Avon. On the other hand 
there is a striking resemblance between these 
measurements and the criterion derived from sub- 
jective impressions and pulse testing in motion 
picture theatres and set up recently by Morr [7] 
at a meeting of the Society of Motion Picture and 
Television Engineers. He claims that echoes more 
than 100 ms after the direct sound should have 
an intensity level at least 20 dB below the direct 
sound. This compares extremely well with the 
derived criterion for speech (Fig. 7a); Morr would 
be mainly interested in speech because as A.W. 
COLLEDGE said at the same meeting “‘about 80 or 
90 per cent of all sound coming out of Hollywood 
is speech”’. It appears likely to the present authors 
that the differing requirements of a motion pic- 
ture theatre and a live theatre lead to different 
design criteria. With motion picture theatres 
where the amplifier gain can be increased to pro- 
vide an adequate listening level disturbing echoes 
need not be tolerated. In a live theatre the audi- 
ence demands an adequate listening level even 
although this may entail the use of definite echoes. 

Finally it should be pointed out that it is fairly 
certain that a still lower reverberation time in the 
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listening room would not change the acceptability 
of the echoes since, assuming the sound to decay 
exponentially, the echoes from the room will be 
far below the tolerable limit for artificial echoes. 
Experiments out of doors by the authors gave 
results similar to those reported here. On the basis 
of results reported above there is a possibility 
that echoes of —35 and —40 dB are tolerable 
at any delay merely because they approach the 
noise level; however actual listening showed that 
inputs of —50 dB were not completely masked 
by noise and hence it must be assumed that the 
—35 and — 40 dB results shown in Fig. 5 can be 
accepted as substantially correct. 


6. Conclusions 


Speech with a single artificial echo presented 
in a non-reverberant room disturbed listeners at 
intensity levels well below those reported by 
Haas. This is due to the lack of room echoes here 
whereas they were present in earlier experiments. 
Louder echoes are acceptable with music than 
with speech. An acceptable-echo-level/delay cri- 
terion deduced from the results is in good agree- 
ment with a criterion developed from pulse test 
work in motion picture theatres where numerous 
echoes to obtain satisfactory listening intensity 
are not required. Similar criteria derived for string 
and organ music allow louder and longer delayed 
echoes and this is in agreement with previous 
criteria of longer reverberation times for music. 
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THE EFFECT OF WALL SHAPE 
ON THE DECAY OF SOUND IN AN ENCLOSURE 
by J.W. HEAD 


British Broadcasting Corporation, Research Department, London 


Summary 


A two-dimensional room is considered in which three of the “walls” are straight and the fourth 
has small, regularly spaced symmetrical projections. The method of FrEsHsBAcu [1] is applied 
to determine the effect of the projections on the natural frequencies and normal modes of the 
room; explicit formulae are obtained for rectangular, circular and triangular projections. 

Recent experimental results obtained by the B.B.C. Research Department [2] suggest that 
studios or halls with rectangular coffering may have superior subjective properties to those having 
other types of projection with comparable dimensions. The results here given appear to provide 
a theoretical reason for this superiority of rectangular coffering. 


Résumé 


On considére une salle 4 deux dimensions dont trois des «parois» sont uniformes tandis que la 
quatriéme porte de petites saillies symétriques 4 espacements réguliers. 
On utilise la méthode de FEsusacu [1] pour déterminer l’effet des saillies sur les fréquences 


propres et les modes normaux de la salle. On donne des formules explicites ayant trait aux saillies 
en forme de rectangle, de cercle et de triangle. 

Des expériences récentes de la B.B.C. Research Department [2] semblent montrer que les 
salles et studios avec saillies rectangulaires sont meilleures 4 un point de vue subjectif que celles 
ou les saillies ont d’autres formes, 4 dimensions comparables. 

La présente étude semble donner une raison théorique pour l’avantage des saillies rectangulaires. 


Zusammenfassung 


Es wird ein zweidimensionaler Raum behandelt, in dem drei der ,,Wande“ glatt sind, wihrend 
die vierte Wand kleine, regelmafig angeordnete symmetrische Vorspriinge besitzt. Mit Hilfe der 
Methode von Frsupacu [1] wird der EinfluB dieser Vorspriinge auf die Eigenfrequenzen und 
Eigenfunktionen des Raumes bestimmt. Formeln werden fiir rechteckige, runde und dreieckige 
Korper abgeleitet. 

Die Ergebnisse von Versuchen, die neulich vom B.B.C. Research Department ausgefiihrt 
worden sind, weisen darauf hin, dafs Aufnahmerdume oder -siile, deren Winde mit rechteckigen 
Platten ausgekleidet sind, nach subjektiver Wertung den Raumen akustisch iiberlegen sein 
diirften, in denen Kérper mit vergleichbaren Abmessungen, aber anderer Form, verwendet 
werden. Die hier angefiihrten Ergebnisse scheinen die Uberlegenheit rechteckiger Platten theo- 


retisch zu begriinden. 


1. Introduction 


The main problem of practical acoustics can be 
briefly stated—to find a way of designing or ad- 
apting studios or concert halls so that they shall 
be regarded as “good” by speakers, musicians or 
listeners who have to use them. Although such 
people may have an immediate subjective ap- 
preciation of the qualities of a particular studio 
or concert hall, it is very difficult satisfactorily to 
correlate this with any objective quantity under 
the control of the designer. Recent experimental 
results obtained by the British Broadcasting 
Corporation Research Department [2] suggest 
that rectangular coffering may be more efficient 
as an agent for improving the subjective quality 
of studios and halls than other departures from a 
smooth wall shape having comparable dimensions. 


As circular cylinders and spheres are widely used 
in the hope of improving quality, it was consid- 
ered advisable to investigate whether there was 
any theoretical reason for this apparent superior- 
ity of rectangular coffering. 

The mathematical formulation of the problem 
in the general case in a practically useful way is 
difficult. We have therefore confined our atten- 
tion to the case of a nearly-rectangular two- 
dimensional room in which three of the walls are 
straight and the fourth has (n+ 1) small sym- 
metrical projections on base d and having a 
maximum depth h, as in Fig. 1. Each projection 
is separated from its neighbour by a straight 
segment of length d, so that the total length a of 
the wall with projections is (2n+1)d; Fig. 1 is 
drawn with n= 3, h,, d are regarded as small 
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compared with a or b, but may have any relation 
to the wavelength. Frsuspacu [1] derives for- 
mulae for the eigenfrequencies and normal modes 


associated with a boundary having irregularities; 


we have considered only the case in which the 
eigenfrequencies are in general distinct. These 


ad td d fd 1 dd td 


A 
ct 


beck a =! 
Fig. 1. General shape of room considered. 


formulae are briefly explained and applied to our 
problem in Section 2. The formulae are perfectly 
general and could theoretically be applied to any 
shape of projection and to unequal projections 
unequally spaced. We have considered in detail 
only three types of projection, namely rectangles, 
circles and triangles; in each case we have taken 
the maximum depth h, to have the same value 
d/2. Explicit results are obtained in these three 
cases; they are given in Section 3. Results are 
discussed in Section 4, their subjective signifi- 
cance in Section 5, and conclusions summarized 
in Section 6. The Appendix, relevant to Section 5, 
discusses the general nature of the decay curves 
associated with one,-two, or several significant 
eigentones. 


2. Fesupacu’s formulae 
The differential equation to be satisfied by the 


sound potential y from which pressure and velo- 
city are derived is 


V*p = (1/c*) (p/ at") (1) 
with the boundary condition 
ep /en=0 (2) 


at all points of the boundary, n being the out- 
ward normal. We thus assume that all boundary 
surfaces are hard. 

When there are no projections, there are a 
number of eigenfrequencies w/2z satisfying 


ep? [e? = (p® x® a”) + (q? 7? /b*) (3) 


where ¢ is the velocity of sound, and p, q are 
integers. We regard these eigenfrequencies as in 
general distinct. With each eigenfrequency is 
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associated a potential 

y= Acos(px2x/a)cos(qzy/b) (4) 
if the axes are taken along two of the straight 
walls and the origin at their intersection. 

We shall denote by @;/22 and w,/2a two 
typical eigenfrequencies; the corresponding in- 
tegers in (3) will be called px, q, and pr, q: respect- 
ively and the corresponding ¢ and A in (4) will be 
called q,, Ax and q:, A: respectively. If 


A= 4/ab 
A2=2/ab 


(p, q both different from zero) 

(5) 
(one of p, q zero) 
the potentials g form an orthonormal set satisfy- 
ing (1) in the area of the room without projections. 
This condition of orthonormality is required for 
the derivation of FEsHBacu’s formulae (6), (7), 
(8) and (9) below. 

Let R denote the region inside the room with 
projections, and let S denote the boundary ob- 
tained by going round the room anticlockwise 
and following the curved surface of the projec- 
tions. Then we define 


Au =o (@p./én) dS, Nu=[grdR. — (6) (7) 


For the room with projections, the eigenfrequency 
Q,/2a corresponding to w,/2z is given by FEsH- 
BACH [1] as 


QF =a + (2 An|Nu) + ct Y Aig|(oi — os) (8) 
k 


the symbol &’ meaning that the term with k=t 
is omitted. The associated potential function is 


Di = Prt PY Apelor ox) (9) 
k 


where ¥Y, equals g, inside S and changes dis- 
continuously to zero on S. 

Hence the extent to which 2,, D, differ from @,, 
YW, (which are the corresponding quantities when 
there are no projections) depends upon the co- 
efficients Aj:. 

We shall see that in the cases for which we 
obtain Ax, explicitly, it is an oscillating function 
of k for a given t; the highest values of | Ax,| 
obtainable for the various types of projection ap- 
pear to be of the same order of magnitude, but in 
determining their effect on (8) and (9) we have 
to consider not only the magnitude of | A,,| but 
also of the denominator (w; — «;) likely to be 
associated with it. An individual A,, may acquire 
unusual importance if associated with an ab- 
normally small value of (@; — w;). We shall con- 
fine our attention to the broad, general behaviour 
of three different types of projection for which 
(6) and (7) can be obtained explicitly, and take 
little notice of isolated abnormal cases. 
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3. Explicit results for rectangular, circular 
and triangular projections 


For any type of projection, we obtain Aj, An 
from (6) and N;, from (7). If 6 is the angle bet- 
ween the outward normal to the projection and 
the negative x-axis, we thus obtain 


Aut = — ¥\ At Ax{ cos (pex/a) cos (query /b) - 
r=0 cS 


-[(pex/a) cos 6 sin (pi x/a) cos (qiay/b) (10) 

+ (q:2/b) sin 6 cos (prax/a) sin (qiy/b)] dS 
where dS is an element of the circumference C, 
of the (r+ 1)th projection, directed so that S in- 


creases with x at all points where the tangent to 
C, is not parallel to the y-axis. This reduces to 


Axt = (At Ax pi/4a)- 
*[p (pit pes e+ Qe) + P( Pst Pes Ge — Qe) + 


+ 9 (Pt— Pe Qe+ qe) + P(Pe— Pes qe — qu) ] + 
+ (A, Axquz/4b) - (11 
‘[p (pet pes qet qx) + p(pe+ Pes qe— qk) + 
+ p(Pe— pes e+ qe) + p(pe— Pes qe — 4) ] 
where 
p(x, 285] { sin (x2x/a) cos (Bry/b) cos 0S, 
r=0C, (2) 


y (a, B) = a J cos («zx /a) sin (Bay /b) sin dS; 


id 


x, B being integers. 
The corresponding result for Ni, is 


Nu=Ai foot ee as | cos? (quay) dy — 
(13a) 
(2r+1)d b 


Ar { cos? (pix/a) dx { cos? (qizy /b) dy 
2rd rs 


ina: °. 


where g—b—f(x—2rd), and f(x—2rd) is the 
depth of the projection at a distance d along it. 
The first integral reduces to 1 in all cases. In the 
(r+ 1)th term of the second we put x=(2r+o)d 
and obtain 

1 Ard: 
4 (13b) 


“J Af (od) + (b/2qu21) sin [(2qu-"/b) f (0€)]} - 


Nu=1l— 


n 


> cos [(27 + 0) (2pixd/a)]} do 


r=0 


-{(n+1)+ 


=1 —A? 7(p., 41), say, (14) 
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where f(x —2rd)=f(cd) is the depth of the pro 
jection at a distance d along it. 

If f(cd) is continuous and differentiable in the 
whole range 0<o <1, and f’ (cd) is always finite, 
then it follows that 


aBy(x,B)=bay(a,f). 


For since § is the angle between the outward 
normal to the projection and the negative direc- 
tion of the axis of x, 


cos# dS=—dy, sin§ dS=dx. 


If in the (r+ 1)th terms of p(a,f) and p(x,f) we 
put x= (2r+o)d, the former becomes an integral 
over the range 0 to 1 with respect to o which can 
be integrated by parts. The integrated term 
vanishes at both limits since f(0)=f(d)=0; the 
residual integral reduces to the (r+ 1)th term of 
y(x,B) multiplied by (bx/af); this proves (15). 
The proof can be extended to the case in which 
f' (cd) has a finite number of isolated infinities in 
the range 0<o<1 provided that f(od) remains 
continuous and differentiable. This extension 
covers the case of circular projections, but not 
that of rectangular projections, where f(ad) is 
discontinuous for c= 0 or 1. 

We now have to apply (12), (13) and (14) to the 
cases of rectangular, circular and triangular pro- 
jections respectively, with due consideration for 
modifications required when (i) p:, Px, qe OF qx is 
zero, (ii) Pk=pz, and (iii) qx= qr. 

In case (i), (11) will be affected because the 
product A,A, will not have its normal value 4/ab 
(see (5)). (12) will not be directly affected unless 
case (i) is combined with case (ii) or case (iii), 
which means that p:=pr=0 or qr=qe=0. Cer- 
tain terms in the normal evaluation of (13b) will 
become indeterminate in case (i) and these should 
be replaced by their limiting values if p, or q is 
regarded as tending continuously to zero. Like- 
wise, in cases (ii) and (iii), certain terms in the 
normal evaluation of (12) become indeterminate. 
The correct value for p(«,f) and w(x,8) can al- 
ways be obtained in such cases by taking the 
limit when « or f is regarded as tending continu- 
ously to zero. We therefore do not need to con- 
sider these special cases further; they do not 
change the general nature of the results, and 
explicit formulae will be given only for p:, px, 
Qt> Qk, P1— px and q:—qx different from zero. 


(15) 


(a) Rectangular projections 
In this case, in (12), 
6=0 dS=-—dy forb=y=b— bd, x=2rd 


6=2/2 dS=dx for2rd<x<(2r+1)d,y=b— ed 
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6—x dS=dy for b—S<y<b,x=(2r+l)d 


and in (13b), 
f(x—2rd)= ‘4d. 


Hence 
plop) =>) J sin (axrdja)cos (Bry /b) dy 
r=0 b—%4d (16) 
= (— 1)’ (b/Bx) sin (Bd /2b) - F (x) 
where 
F(a) 2 x sin (xard/a) = a” sin [«ar/(2n+ 1)] 
=cot[«7/2(2n+1)] (a odd) (17) 
or 0 (x even). 
Also 
n (2r+1)d 
wp (c,B) = fies i cos («2x /a) - 
-(— 1)’ sin (Bd /2b) - dx (18) 
= (—1)°** (2a/xz) G («) sin (xd /2a) - 
‘sin (B2d/26) 
where 
G(x) = x cos [(2r + %) aa/(2n+1)] 
= 0 (« odd) (19) 


or % sec [az /2 (2n+1)] 

(x even, not a multiple of (2n+1)) 
or (n+1) cos [xz/2(2n+1)] 

(x a multiple of 4n+2). 


Finally, 


1 (Pog) = (4/4) H [(d/2) + (b/2quzx) sin (q.zvd/b)]- 


. +) + eos [2 piad (2r +) /al} do 
= (1/8) [d + (b/q:) sin (q:d/b)]- (20) 
*[(n +1) d+ (a/pix) sin (pi ad /a) G (2p:)]. 


(b) Circular projections 
In this case, in (12) we put 
x=2rd+ %d(1—cos6), 
y=b— \dsin§, 
dS= %d- dé 


(21) 


and the resulting integrals can be evaluated with 


the help of the identities [3] 
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en, #=J,(2) +2 >, Jz,(z)cos2p0 + 
p=l 
4+ 2i >) Jep—1(2) sin (2p’ — 1) 6, 
Pit 
J,+1(2) + J,—1(2) = 27 J, (2) /z, 
(22) 


Joe 2S n21(z), = — 25, (2) 


where the J’s are Bessel functions of the first kind, 
and » is arbitrary. 

Although sin(azx/a), cos(asx/a), sin(By/b) 
and cos(fzy/b) thus give infinite series in sines 
and cosines of multiples of 6, only a few terms 
contribute to the integrals (12) transformed by 
(21). We find 


aBy(x,B) =bay(«,p) =ab(—1)'H(x,f) (23) 
where H(«,f) = 2G (x) s (ele (2p 1)s 
ee (24) 


*Jop—1 («2d /2a) Jop_1 (Bad /2b). 
To evaluate y(p:,q:) in (14), we put 
o=%(1—cos6), f(od)= %dsin# 


and obtain 


mt 


z(Psa)= Ig | [dsind + la) sin (q:2dsin6/b)]- 


. [on 1)+ a cos [(4r + 1) pixd/a]- 


-c08(prdeos6a)| sin 6 dé 
a= e [ % (n+ 1) ad + (n +1) (b/q:) J, (qd /b) + 
+ (a/ps) J, (perd/a) G (2p.)] + 
+ (ab [8c peqs) H (2p, 241) (25) 
(c) Triangular projections 


In this case, in (12), for 2rd<x<(2r+ %)d, 
we write 


6=n/4, x=2rd+ %éd, y=b— %éd, 
dS—d)/2dé/2, 
while for (2r+ %)d<x<(2r-+ 1)d, we write 
6=3n/4, x=(2r+1)d—%éd, 
y=b—%éd, dS=d/2dé/2. 
We find 
(8) =(-~ 1)" (x) abtaja(arp®—be4])- 
- [cos (x2d/2a) — cos (B2d/2b)] 
and (15) applies, so that 
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p(x, B) = (aB/bx) p(x, B) (27) 


while to evaluate 7(p:,q:) we have 
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Table I 
Maximum values of p(x,/) and p(«,f) 


/ Type of Maximum] Value of x:8 |Maximum| Value of «:f 
= SoG ’ E 
fled) oa 0S05% projection | p(x, B)| ly(a,B)| | for maximum 
f(od)=(1-o)d ¥%soel 
i z 2n+ 2); 
so that Rectangular 1; Bd/b<a/2| 2a/z? eer) ) 
, = 2) | 14 (n +1) d2 x? 28 : 2n or 2n+2; ; 2n or (2n+ 2); 
4(Pe qe) = (1/82 { % (n+1)d?2?+ (28)  Girenlar pd/b << 7/2 0.66a/7 | ear 1.18b/d 
+ (n+ 1) (b?/q,)[1 — cos (q:2 d/b)]+ ef onl be dpateetic le Qn orale 
ee Bd/b-<x/2 = near b/d 


+ (a°/p?) G (2p.) [1 — cos (pixd/a]+ 
+028? G (2p,)[(a?q — 6*pi) 
- [cos (piad/a) — cos (qird/b)]}. 
4. Discussion of results 


The outstanding features of these formulae for 
the various types of projection are: 


(i) For circular and triangular projections, but 
not for rectangular projections 


aB yp (x, B) = bay (x, f), (29) 
Ax. = 0 (30) 
if one of p:, px is odd and the other is even. 
(ii) For rectangular projections 
¢(x,p)=0 (31) 
if p:, px are both odd or both even, 
y(a,)=0 (31) 


if one of p:, px is odd and the other is even. 


(iii) If the wall without projections has any 
abnormal concentration of eigenfrequencies, the 
projections will tend to break up this con- 
centration. For if in (8), gy, and gx are different 
but wx, the term c44j,/(@;—«;) in (8) will 
be unusually large and 2, will therefore differ 
appreciably from w,;. Similarly 2, will differ from 
cw and therefore from @;, and these differences 
will not in general be equal. From (9) we see that 
inside the room, ®, will be mainly a combination 
of gy; and g:, while ® will be mainly a different 
combination. We have here assumed that some 
of the Ax; associated with the abnormal con- 
centration of eigenfrequencies are sufficiently 
large. A set of projections such that all the Ag, 
associated with a particular concentration of 
eigenfrequencies were zero or very small would 
fail to break up that particular concentration. 


(iv) The maximum numerical values of g(x,/) 
and y(«,f) for the various types of projection, 
and the values of x,f for which they occur, are 
very approximately calculated and tabulated 
below: 


It is clear from (30) and (31) that if one of ps, px 
is even and the other odd, only rectangular pro- 
jections can give a non-zero value for Ax. The 
essential difference for our present purpose bet- 
ween rectangular projections and other types is 
that rectangular projections have finite portions 
which are perpendicular to the wall on which they 
are based, whereas circular projections have only 
infinitesimal portions perpendicular to this wall 
and triangular projections have no such portions. 

The maximum values of g(«,f) and p(x.) are 
of the same order of magnitude for the three 
types of projection, but small denominators 
(w;—«;) are likely to arise when p:~p, and 
qi qx so that p(p:—px: q:—qe) would make a 
relatively large contribution to A,, for rectangular 
projections if (p:—ps) is odd and small, and 
(q: — qx) is small. Small values of (@; — w;) ean also 
arise in association with even values of (p;+ px) 
and (q:+ qx) which would make |g(x,f)| or 
|w(x,$)| near a maximum in the other cases, but 
this is much less likely. For non-rectangular pro- 
jections, small (; — @;) will normally be associat- 
ed with very small Ax:. 

It thus appears that in the case of rectangular 
projections the terms Ax:/(w; —@;,) will make a 
more significant contribution to (9) than in either 
of the other cases. 

If the projections were not evenly spaced, or 
occurred on more than one wall, or the room was 
three-dimensional, the above procedure could in 
theory be applied, but would in practice be im- 
mensely more complicated. We should not then 
expect half the A,, to be zero for circular and 
triangular projections as in the case we have con- 
sidered, but we can reasonably expect that the 
contribution of the portions of rectangular pro- 
jections which are perpendicular to the wall from 
which they project is to increase sharply many of 
the Ax, which would be small for the other types 
of projection, and especially to increase certain Ag; 
associated with small (w;—w;). We are thus led 
to expect more significant terms in (9) from rect- 
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angular projections than from other projections 
similarly distributed and of comparable dimen- 
sions. 


5. The subjective significance of the results 
The B.B.C. Research Department [4], [5], [6] 


have found that serious faults in studios and halls 
are often associated with recognisable character- 
istics of the pulsed-glide decay curves for a 
limited range of frequencies. In these decay 
curves, the logarithm of the envelope sound pres- 
sure during decay is plotted against time. For a 
subjectively good studio or hall, the successive 
pulsed-glide curves obtained as the frequency of 
the source is slowly increased (but the positions 
of source and microphone are fixed) should not 
generally show marked irregularity (departure 
from linearity) or marked lack of irregularity, and 
they should not form clearly-defined patterns. 

We consider in the Appendix the theoretical 
decay-curve which would be obtained on the as- 
sumption that all but n of the eigentones excited 
by the source are negligible. We there show that 
the irregularity of the decay curve is connected 
with the number and relative importance of the 
significant eigentones; decay curves of moderate 
irregularity occur when the two principal eigen- 
tones are comparable without being of nearly 
equal intensity, or when there are several signi- 
ficant eigentones of comparable intensities. Equa- 
tion (9) shows that the effect of any type of pro- 
jection is to replace the potential function gy; 
(associated with a particular eigentone of the 
enclosure without projections) by a linear com- 
bination of neighbouring eigentone potentials, 
and therefore in general to reduce the relative 
importance of eigentones which were outstanding 
for the enclosure without projections. Projections 
should therefore improve the subjective quality, 
and rectangular projections (for which, as we 
have seen, the difference between ®; and Y; in (9) 
is most marked) should improve the subjective 
quality more than other types. 


6. Conclusions 


(a) Rectangular projections alter the natural 
frequencies and associated normal modes of a 
smooth-walled rectangular hall or studio more 
than circular or triangular projections having the 
same base and maximum depth. 

(b) This result is due to the fact that rectan- 
gular projections, unlike circular and triangular 
projections, have finite portions perpendicular to 
the wall on which they are based. 

(c) Rectangular projections are also more ef- 
fective for breaking up abnormal concentrations 
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of eigenfrequencies than circular or triangular 
projections. 


(d) We have only considered the effect of pro- 
jections on normal modes in general; for a-few 
isolated particular normal modes the rectangular 
projections might happen to be relatively in- 
effective—for example, if q: and q, are both mul- 
tiples of 2b/d. A slight alteration of d or b would 
be adequate to correct this ineffectiveness in the 
rare cases when correction was necessary. Each 
type of projection will have a different series of 
modes for which it is thus specially ineffective. 


(e) The general effect of projections would ap- 
pear to be the reduction of the prominence of the 
leading eigentones at most frequencies. This effect 
would be most marked where one eigentone is 
outstanding, or where only two eigentones of 
nearly equal intensity are significant; no marked 
effect would be expected from the projections 
when several eigentones are of comparable im- 
portance. Serious faults in actual studios are 
often correlated with the type of pulsed-glide 
curve associated with the predominance of one 
eigentone, or of two eigentones having nearly 
equal intensity, so that any form of projection 
would be expected to improve the subjective 
effect at most frequencies. 


(f) The general effect of projections being an 
overall subjective improvement (conclusion (e)) 
and the effect of rectangular projections being 
more marked than that of other types of compar- 
able dimensions (conclusion (a)), we should expect 
the overall subjective effect of rectangular pro- 
jections to be superior to that of projections of 
other types. 


(g) Although the case here considered mathe- 
matically is necessarily a special one, it should be 
adequate as a qualitative first approximation for 
other cases in which we should expect similar 
results but the mathematical formulation of those 
results would be very difficult. The theory here 
set out is to be regarded as a first approximation 
rather than as a final result, and is put forward 
tentatively in the hope that it may point the way 
to a better approximation. 
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Appendix 
Theoretical pulsed-glide decay curves 


when n eigentones are significant 


We suppose that when the source ceases at 
t=0, n eigentones are significantly excited. Let 
the initial amplitude, damping rate, frequency 
and phase of the ith eigentone be respectively 
Ci, wi, wi/2a and qi. Then the resultant sound 
pressure P at time t is given by 


J ics 2 Ciexp[— wit+i(wit+qi)]. (1) 
We assume that the C; are in descending order of 
magnitude, and that all differences w:i—«j; are 
small compared to the w; themselves. The quan- 
tity measured by the pulsed-glide technique is 
proportional to log R, R being the magnitude of 
the vector P determined by (1), so that 


R=) Cee +2) 2 GG: (2) 
i= i=l j=i+ 


“exp [— (ui + mj) t] cos [(i — @j) t+ (pi — 9y))- 
If n=1, R=C,-e “', so that the pulsed glide 


is a straight line. 

If n=2, and py=p, Re oscillates with 
period 22/(w,—@,) [s] between C,—C, and 
C,+ C,. Hence if C,=C,, the pulsed-glide will 
have a dip to —oo followed by a maximum in 
each beat cycle interval of 2/(w,—@,) [s]. If C, 
is less than, but nearly equal to, C,, the pulsed- 
glide will have a maximum and a minimum in 
each beat cycle, whereas for smaller, but not 
negligible, values of C,, the pulsed-glide will mere- 
ly have an inflexion in each beat cycle. If “4, < v5, 
the pulsed-glide will for sufficiently large C,< C, 
have maxima and minima in each beat cycle at 
first, then mere inflexions, and finally the pulsed- 
glide will become straight with slope correspond- 
ing to the damping rate jy. If > uw, Cy> Cy, we 
can have a pulsed-glide which starts by being 
straight with slope corresponding to damping 
rate “4, then after a few beat cycles inflexions ap- 
pear, then maxima and minima, then inflexions, 
and finally a straight portion of slope correspond- 
ing to damping rate jw, (the so-called “double- 
decay” type of pulsed-glide associated with sub- 
jectively undesirable effects). If the relative phase 
of the two eigentones happens to be right, R may 
pass through a single zero during the time of 
maxima and minima. Some of the above features 
may be lost in practical pulsed-glides because R 
falls below the detectable threshold before they 
have developed. 

We can sum up by saying that if n= 1, i.e., one 
eigentone is outstanding, the corresponding 
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pulsed-glide is very straight. If two eigentones 
predominate above all others and are unequally 
damped in such a way that they are of equal 
intensity at a particular instant T during the 
decay, the corresponding pulsed-glide is very ir- 
regular near time T. The extreme case is that of 
two eigentones of equal intensity and equally 
damped, for in this case there is in effect a time T 
once every beat cycle. 

If several eigentones_are significant, it is still 
possible, but unlikely, for R to be zero. The con- 


ditions for R to be zero are that simultaneously 


R= DY Cie“ cos[(wi — ws) t+ (vi—;)] =9, (3) 


=! 


i 2d, Ci-e “i sin [(w:—wj) t+ (pi—qy)] =0. (4) 


There will at most be a finite set S, of values of t 
for which (3) is true before R falls below the 
threshold. Similarly there will be at most a finite 
set S, of values of t for which (4) is true before R 
falls below the threshold. In general the sets S, 
and S, will have no common member. Either or 
both may have no members. Thus in general we 
should not expect R to have a zero when several 
eigentones are significant though it is possible in 
particular cases for R to have either a zero or a 
small value. Extreme values of R require the 
satisfaction of stringent conditions by the relative 
phases (wi—wj)t+ (pi—q;). If we regarded these 
relative phases in (2) as completely random, we 
should expect R to be near 


Ry = | pa Cie, (5) 


In general, therefore, when there are several 
significant eigentones of comparable intensities, 
we should expect R to oscillate moderately about 
the value Ry, so that a pulsed-glide of moderate 
irregularity would be obtained. 


References 


[1] Fesupacu, H., On the perturbation of boundary con- 
ditions. Phys. Rev. 65 [1944], 307—318. 

[2] SomerRvILLE, T. and Warp, F.L., Investigation of sound 
diffusion in rooms by means of a model. Acustica [1951], 
40—48. 

[3] Ancor, A., Compléments de mathématiques. C, N. E.T., 
Paris 1949, pp. 361, 368. 

[4] SomervittE, T., Acoustics in broadcasting. Rep. of 
Building Research Congress, London 1951, pp. 53—59. 

[5] Somervitte, T. and Grtrorp, C.L.S., Composite 
cathode ray oscillograph displays of acoustic phenomena 
and their interpretation. B. B.C. Quart. 7 [1952], 1—13. 
This is substantially the same as: 

[6] SomerviL_e, T. and Grirorp, C.L.S., B.B.C. Pulsed- 
glide acoustic displays, parts 1—3. FM-TV Radio Com- 
mun. 12 [1952], No. 6, 22, 43; No. 7, 28, 30; No. 8, 22, 38. 


REMARKS ON THE CONCEPT OF ACOUSTIC ENERGY 


by A. SCHOCH 


Institut fiir theoretische Physik der Universitit Heidelberg 


Summary 


The conventional expressions for acoustic energy—recently found not to be correct in a strict 
sense—are shown to have characteristic properties by which their use in acoustics can be justified 


satisfactorily. 
Sommaire 

On montre que les expressions usuelles pour l’énergie acoustique — récemment prouvées 
incorrectes d’un point de vue rigoureux — ont des propriétés caractéristiques, qui justifient 


parfaitement l’usage de ces expressions dans l’acoustique. 


Zusammenfassung 


Die tiblichen Energieausdriicke der Akustik — von denen in letzter Zeit gezeigt worden ist, 
daB sie nicht streng richtig sind — besitzen charakteristische Eigenschaften, durch die sich ihr 
Gebrauch in der Akustik befriedigend rechtfertigen laBt. 


1. Introduction 


The expressions in common use for the energy 
density EF and the energy flux density J of a sound 
wave are 


E=5 & tk (p Pa) = 
Me a. (P 70) ad 
Ae E * (ee)? 
T=(p—p)u, (1.2) 


where u= particle velocity, p—pressure, 0 = den- 
sity, «= 1/oc?= adiabatic compressibility of the 
medium, c= velocity of sound; the suffix 9 indi- 
cates values at equilibrium (medium at rest). The 
expressions are of second degree in the field vari- 
ables u, p —po, and usually correctness to second 
order is assumed, if solutions of the first order 
approximation to the acoustical equations are 
substituted for u and p— po. 

The expressions (1.1), (1.2) have been sub- 
jected to criticism by ANDREJEw [1] and recently 
by J. J. Markuam [2], the objections being that 

(i) in most of the textbooks it is not specified, 
whether the expressions refer to Eulerian or 
Lagrangian representation of the motion, i.e. 
if E and I refer to volumes and areas fixed in 
space, or to volumes and areas consisting al- 
ways of the same particles and accordingly 
changing during the motion; 

(ii) second order terms have been omitted in E 
and I. Unfortunately, to take account of the 


missing terms would require solutions of the 
acoustical equations to the second order of 
approximation. 


The present note is an attempt at a rehabili- 
tation of the expressions (1.1), (1.2) by showing 
that the continuity equation relating the total 
energy and the total energy flux can be separated 
into two independent equations, one of which 
relates the quantities of (1.1), (1.2)}. The latter 
of these equations is of course well known, though 
doubts had arisen as to its interpretation [2]. 


2. Energy equations of hydrodynamics 


The following considerations refer to a liquid 
(or gaseous) medium without any dissipative 
properties (as for instance viscosity, thermal con- 
ductivity, structural relaxation phenomena, etc.). 
In such a medium the state of stress is given by a 
scalar pressure p, which is furthermore complete- 
ly determined by the density 0 (resp. the volume 
1/o per unit mass) of the medium: 


P= P(e). (2.1) 


The hydrodynamical equations of motion, to- 
gether with the equation expressing conservation 
of mass, immediately lead to the following 
equation: 


hg (‘“) ae (*)| + (pu) =0. (2.2) 


1 This separation has already been used in a paper on 
the momentum of a sound wave [3]. 
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Vv refers to fixed coordinates and 


= = x +(uV) (2.3) 
gives, in the usual notation, the rate of change of 
a quantity attached to a particle, when following 
that particle in motion. 
pu= I: the flow of energy through a material 
unit surface (i.e. p(uda), is the work done per 
unit time on some part of the medium by an ad- 
jacent part in contact along the surface element 
da). u2/2 is the kinetic energy per unit mass, and 
Dil DU 
tsi] he 
is to be interpreted as the rate of change of an 
internal, potential energy U per unit mass of the 
medium. The existence of a potential energy 
follows from the assumed relation (2.1) between p 
and 0: 


l/o 


lal) 


1/05 


U is the work done by compression of the unit 
mass from volume 1/o to 1/0). Apparently the 
reference value 0), and with it an additive con- 
stant in U, are arbitrary. 

By introducing U, (2.2) takes the form 


D (uw Boon 
0D, (= arpr) sp Fesey, (2.4) 
Applying (2.3) and using the conservation of 

mass leads to the equivalent equation [1], [4] 
ck 


= + V (1+ Eu) = 0, 


(2.5) 


with 
2 


E=o (5 ae u) = energy density. (2.6) 


The meaning of these equations is conspicuous: 
The balance of energy is described by (2.5) for an 
element of volume fixed in space, i.e. in the 
Eulerian sense (in this case the energy flux is 
augmented by the convective energy transport 
Eu), and by (2.4) for a material element of the 
medium, consisting always of the same particles, 
i.e. in the Lagrangian sense. (The \/-operations 
occurring have to be performed in an Eulerian 
representation of the field variables in both cases.) 

The energy theorem (2.4) can, in its substance, 
be found in Lams’s Hydrodynamics [5]. A juxta- 
position of the alternative forms (2.4), (2.5) of it 
seems to have been given for the first time by 
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ANDREJEW. The equations have been repeated 
here for convenient reference. 

As was mentioned above, an additive constant 
is arbitrary in U. In E this corresponds to an 
arbitrary additive quantity, which is propor- 
tional to 9. For this reason, it is not possible to 
fix the scale of energy density by taking the 
difference of the energy density in the state in 
question and in some reference state (as ANDRE- 
JEW [1] did in his definition of Eulerian energy 
density). In fact, the conservation law (2.5) would 
be violated by adding a constant to E. Physically, 
this means, that potential energy is attached to 
matter, like internal energy in thermodynamics 
and like kinetic energy. Things may, however, 
become more complicated if longe range interac- 
tion forces play a part, as e.g. in the presence of 
electric polarization (piezoelectric media). 


3. Decomposition of hydrodynamic energy 


The potential energy per unit mass can be 
written 


rj 


1/og 


1 ite 

me (- = (Pp ~ Po) d (") = ta 
i S 1/05 a 
Ur Un 

= Uy+ Un, 


where py the pressure in the state of rest. Then 
Uir is, to second order terms, equal to the con- 
ventional acoustic energy density (1.1). U1 is in 
many texts supposed to vanish in the time aver- 
age. Usually, however, a non-vanishing second 
order term of Ut is left even in the time average, 
as has been shown for the case of a plane pro- 
gressive wave [1], [2]. Therefore the energy of a 
sound wave in a strict sense comprises both 
Ur and Un. 
Carrying out the corresponding separation 


aly ea Aaa, 2 (; 


0) — (P— Po) 4 (;) 
in equation (2.2), we derive 
epe(s) Poel) Pa pele) # 
+ V (pou) + V[(p — po) u] = 0, (3.2) 


and as 
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on account of the equation of conservation of 
mass, (3.2) separates into two equations 


D Ur = 7 i 
@—, + ee, pi, (3.3) 

D (w - oi = 
onlst Un) +vIn=0, Tii=(p—Po)4 (3-4) 


which could also be put into the equivalent Eule- 
rian forms corresponding to (2.5). 

(3.4) is the conservation theorem usually deriv- 
ed from the first order approximation equations 
of acoustics and provides the basis for the con- 
ventional energy expressions (1.1) and (1.2). With 
the definition (3.1) of Un it holds rigorously, and 
not only within the framework of the linearized 
acoustical equations. 


4. Discussion of Ui and Ik 


As Us, It contain the field variables linearly, 
solutions to the second order of approximation 
are required to give values of the energy quan- 
tities correct to second order. In principle, first 
order solutions would therefore be insufficient for 
the computation of acoustic energy. 

It seems, however, that the contributions Ur, Ti 
are far less important in acoustics than Un, Tu. 
To understand the physical meaning of the 
I-terms, let us compute the total I-energy con- 
tained in a fixed quantity of matter (dr = element 
of volume): 


pie tT 
fetrae=—m [ o(o—2)) dt = — po(V— Vo). 


Here V and V, are respectively the instantaneous 
and the equilibrium volumes of the quantity in 
question. Apparently the total I-energy equals 
the work done by the equilibrium pressure p, 
during the change of the volume. Therefore, the 
total 1-energy is zero, if there is no change of the 
total volume. This condition would be fulfilled 
for a medium contained in a rigid vessel, a case 
which may, however, not be regarded as strictly 
realizable. 

But the following is a case of practical interest : 
Suppose a source of sound waves (e.g. a vibrating 
piezoelectric disc) positioned in an extended 
medium and put into action for a limited time. 
After some time, the sound waves cover a region 
as sketched in Fig. 1. Then the I-energy in a part 
of the medium entirely containing the wave 
packet (i.e. in a part bounded by the surface of 
the source and any material surface 5S as indicated 
in Fig. 1) will be zero, if the source assumes its 


A.SCHOCH: THE CONCEPT OF ACOUSTIC ENERGY 183 


initial volume after having been switched off. 
That is, Ur does not contribute to the total energy 
of the sound waves. 

In a steady-state sound field, maintained by 
periodic vibrations of the surface of the source, 
we similarly find that the energy flux T; does not 
contribute to the energy radiation from the 
source. The I-energy traversing a material sur- 
face S enclosing the source is per unit of time 


fe popes DV 
ppo(udo) = Po aiYa 


if V, for the time being, denotes the instantaneous 
volume enclosed by S. The time average of this 
(in the Lagrangian sense) becomes zero because of 
the assumed periodic motion. The absence of 
radiation due to I-energy might have been in- 
ferred from the case of a finite wave train con- 
sidered in the preceding paragraph: There is no 
net transport of I-energy because the total amount 
of it contained in the wave train is zero. 


Oe tee 


Fig. 1. A wave packet enclosed within a surface S on which 
the medium is at rest. 


The statements may be illustrated by the case 
of a plane progressive wave, where an exact solu- 
tien is available. If the motion is taking place in 
x-direction and is excited by a piston at x= 0, the 
displacement & of which is varying according to 
&= f(t), the solution is most conveniently written 
in Lagrangian representation (i.e. « denoting the 
coordinates of the undisplaced particles). Ap- 
proximated to second order it reads [6], [1], [7]. 


[3]: 


x xXx x\]? 
“i vein gl Ve ic see dl Cd LC 
Het) =F (t=) + 35 Ei ( “a (4.1) 
if the initial conditions are imposed: 
git) —0,. f'(t) = 0 for t= 0 
0& 
E(x,0)=0, —=u(x,0)=0 forx>0. 


ct 
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Herein is c = /dp/dg the instantaneous velocity 

of sound, «=1+ ge a , and the suffix y indicates 
Cy \de), 

values in the equilibrium state. 

With the piston returning to its initial position 
after having excited a sound wave, (4.1) shows 
that also every particle returns to its initial posi- 
tion after the passage of the wave train?. The 
volume of that part of the medium which is 
covered by the wave train is therefore equal to 
the equilibrium volume. Thus the total I-energy 
is zero, as stated above, and as had been shown 
earlier [3]. 

From (4.1) one easily deduces that 


1 Wh hal! 0€ 


The expression in terms of the particle velocity u 
shows that the (Lagrangian) time average of U1 
does not usually vanish, in spite of zero values of 
the time average of u and of the total I-energy. 

That there is—in the time average—no I-en- 
ergy flow in the case of periodic motion is im- 
mediately evident from the vanishing time aver- 


age of pou. 
5. Conclusion 


On account of the properties of Ur and Ir and 
under the conditions prevailing in acoustics, it 
would therefore seem a reasonable convention to 


— 


consider the contributions 0 (< aa Un| = Ey and 

* This would not be true if f(t), f’(t) were not simultane- 
ously zero for t=0. A start of the piston with a finite 
velocity —which would produce a shock wave —shall, how- 
ever, be excluded here. It could not, moreover, be realized 
with a finite mass of the piston. 
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Thr as the acoustic energy density and energy flux 
density proper. Err has the property of never 
being negative, whereas oUr has not. (That 
Un =0 follows from the hydrodynamic stability 


of the medium, which requires ap/a(" | <0.) 

In a dissipative medium, a potential energy no 
longer exists. Equation (2.4) has then to be genera- 
lized by adding the heat flow to J and inter- 
preting U as the internal energy in the sense of 
thermodynamics. But the separation of this gene- 
ralized equation into two equations correspond- 
ing to (3.3) and (3.4) is still possible, with the 
same expressions for U1, Tr as above, satisfying 
the same conservation equation (3.3). It is the 
generalization of (3.4) which contains the thermo- 
dynamics. This means that dissipation (sound 
absorption) feeds from E11, whereas (3.3) de- 
scribes only some reversible shifting about of 
mechanical energy. This consideration may serve 
as further argument for the distinction of part 11 
as acoustic energy. 

The present note was stimulated by a discus- 
sion with Dr. P. J. WESTERVELT, who directed 
the author’s attention to the papers by J. J. 
MarkHAM. The author would like to thank 
Dr. WESTERVELT and Dr. Marxknam for discus- 


sion and correspondence. 
(Received 10th November, 1952.) 
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zur Aufgabe gemacht. Eine Reihe namhafter Fachleute aus 
Wissenschaft und Industrie steht der Gesellschaft als 
korrespondierende Mitglieder mit Rat und Tat zur Seite. 
Neben Problemen der Grundlagenforschung sollen zunichst 
insbesondere solche der Verfahrens-, Verkehrs- und Ferti- 
gungstechnik behandelt werden; mit biologischen und 
medizinischen Forschungsstellen besteht gleichfalls eine 
enge Zusammenarbeit. W. Weber 
27. April 1953 


ZUR AKUSTISCHEN ZAHIGKEITSGRENZSCHICHT 


Von E. MEYER und W. GUTH 


III. Physikalisches Institut der Universitat Gottingen 


Zusammenfassung 


Es wird die akustische Zahigkeitsgrenzschicht untersucht, indem in fortschreitenden Schall- 
wellen die Luftbewegung in der Nahe einer starren Grenzflache durch beleuchtete Schwebe- 


teilchen mikroskopisch sichtbar gemacht wird. 


Summary 


The acoustic viscous boundary layer is investigated for sound waves in air in the neighbourhood 
of a rigid wall using illuminated floating particles seen under a microscope. 


Sommaire 


On a étudié la couche limite visqueuse acoustique en rendant visible au microscope, par 
Véclairage de particules en suspension, le mouvement de l’air, dans une onde sonore progressive, 


au voisinage d’une surface rigide. 


1. Einleitung 

Die Luftteilchen in einer fortschreitenden 
Schallwelle fiihren lineare Schwingungen in der 
Schallausbreitungsrichtung aus. Bringt man eine 
ebene starre Flache parallel der Ausbreitungs- 
richtung in das Schallfeld, so wird in gréBerem 
Abstand von dieser die Bewegung der Luftteil- 
chen unbeeinfluBt bleiben. Ihre Schwingungs- 
amplitude ist unabhangig vom Abstand konstant. 
An der Grenzflache selbst haftet die Luft, und 
infolge ihrer Zahigkeit wird ihre Bewegung in der 
Nahe der Grenzflache gehemmt. Man wird also 
eine der Abb. 1 entsprechende Geschwindigkeits- 
verteilung erwarten, wobei die Lange der Pfeile 
die jeweilige Geschwindigkeitsamplitude andeu- 
tet. Sie ist in groBem Abstand konstant und ver- 
ringert sich zur Grenzflache hin. 

Die mathematische Beschreibung des Vor- 
ganges findet sich bei Kircnyorr, RAYLEIGcH, 
CRANDALL und CREMER [1]---[4]. 

Die aus der Herleitung der Schallwellenglei- 
chung gelaufige Newtonsche Kraftgleichung 


QatopGXP. 
eee ae SOE i Ox 
eu Fee wird durch einen Summan- 
ees den erweitert, der der Za- 
——— higkeit Rechnung tragt 
——__—_— Ov ov e 
I See ane ee eP 
tiie 5 ES 15y2t So Ct ag (2) 


3 =— Grenzflache Darin sind v die Schnelle in 

ht dey der x-Richtung (Abb. 1), p 
Vermutete SaluicHovertel- der Druck, Qo die Dichte, 
Jung iiber derGrenzfliche. 7 die Zahigkeit der Luft. 


Die Vereinfachung in der Behandlung der 
Gl.(1) besteht in der Annahme, da der Druck p 
unabhangig von der Hohe y iiber der Grenzflache 
ist. 0/0t wird durch iw ersetzt. Aus Gl.(1) wird 
dann: 


Thre Lésung lautet: 


LAP f wala 


1@ Oo ox 


aon f - ots"), 


Vv, ist der Betrag der Schnelle im freien Schallfeld. 
Die Schnelle im Abstand y von der Grenzflache 
ist nach Gl. (3) (4) 


v(y) = 


-/#» DO - [sy 
ne Galen Vn° cos ae yte | 2n ’ 


Die zur GI.(4) gehérige Kurve ist in Abb. 2 
skizziert. Sie oszilliert mit einer exponentiell ab- 
klingenden Amplitude um die Gerade |v| = vp. 
Den in Abb. 2 mit A=22(mwo/2n)—*% bezeichneten 
Abstand nennt CREMER die ,, Zahigkeitswellen- 
lange*‘ oder die ,,akustische Grenzschichtdicke*. 
Ersteres, da ja die Gl.(3) formal eine Welle, die 
Zahigkeitswelle, beschreibt, die von der Grenz- 
flache her in der y-Richtung fortschreitet.t Die 
Zahigkeitswellenlange wachst also, wie man aus 
Gl.(3) und im folgenden aus Abb. 6 sieht, um- 


1 Die Analogie zum elektrischen Skineffekt ist offen- 
sichtlich, der ja auch auf einer in das Metall eindringenden 
Welle beruht. 
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Abb. 3. Skizze des Kurven- 
verlaufs zu Gl.(5). 


Abb. 2. Skizze des Kurven- 
verlaufs zu Gl.(4). 


gekehrt proportional der Wurzel aus der Fre- 
quenz. 

Die Phasendifferenz, mit der ein Luftteilchen 
in der Héhe y gegen ein senkrecht dariiber 
im freien Schallfeld befindliches Luftteilchen 
schwingt, ist nach GI.(4) gegeben durch 


Die zu GIl.(5) gehérige Kurve ist in Abb. 3 


skizziert. 


2. Aufbau der Apparatur und Messung 


In einem langen Rohr mit rechteckigem Quer- 
schnitt wird vom linken Ende her die im Rohr 
enthaltene Luft durch einen Lautsprecher zu 
Langsschwingungen angeregt (vgl. Abb. 4). Die 
vom Lautsprecher abgestrahlte Schallenergie 
wird am rechten Ende des Rohres durch einen 
2 m langen Schluckkeil aus Schlackenwolle ab- 
sorbiert, so daB sich zwischen Lautsprecher und 


pape 


/ /~ Bogenlampe 
is oi 
ee 
PUA 
Vv “ Blende 
Facher —\— 
“if 


/ 


77 : 
a8) Grenzflache Absorptionskeil 


Mikroskop 


Lautsprecher 
Abb. 4. Skizze der Apparatur. 
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Keil eine fortschreitende Schallwelle ausbildet. 
Der Rohrquerschnitt (5 x 7 em?) ist so klein, daB 
keine Querwellen entstehen kénnen. Zwischen 
Lautsprecher und Keil sind in den beiden Seiten- 
wanden Fenster angebracht, zwischen denen ein 
Glasréhrchen” von 2 cm Durchmesser liegt. Seine 
auBere Oberflache dient als Grenzflache, tiber der 
die Luftbewegung beobachtet wird. 

Die Schwebeteilchen sind winzige mittels eines 
Zerstaubers in die Luft-geblasene Oltrépfchen. 
Sie werden durch eine Dunkelfeldbeleuchtung im 
Mikroskop sichtbar gemacht. Als Lichtquelle 
dient eine Beck-Bogenlampe, die im Augenblick 
des Photographie- 
rens stark itberla- 
stet wird. Der Spie- 
gel der Bogenlam- 
pe bildet die Licht- 
quelle auf eine 
Blende ab, vor der 
zur Phasenmes- 
sung ein durch 
einen Synchron- 
motor betriebener 
Facher den Licht- 
strahl  periodisch 
kurz unterbricht. 
Ein Beleuchtungs- 
objektivy bildet den 
Facher in die Ob- 
jektebene iiber der 
Grenzflache ab. 
Das von der Licht- 
quelle ausgeleuch- 
tete Gebiet iiber 
der Grenzflache wird durch ein Mikroskop be- 
obachtet. 

Die Oltrépfchen (etwa 10-*°cm Durchmesser) 
schwingen unter dem EinfluB des Schallfeldes hin 
und her und machen die Schallelongation der 
Luftteilchen véllig mit. Ihre Bahnen sind im 
Mikroskop als horizontale Striche sichtbar. Mit 
Hilfe eines elektrischen Feldes (etwa 8000 V/cm) 
kénnen die Trépfchen in vertikaler Richtung auf 
die Grenzflache zubewegt werden, da sie durch 
den ZerstéubungsprozeB fast immer elektrisch 
geladen sind. Man sieht dann im Mikroskop die 
Striche sich auf die Grenzflache zubewegen 
und dabei kiirzer werden, bis sie als Punkt auf ihr 
enden. Macht man von dem ganzen Vorgang eine 
Zeitaufnahme, so kommt die Aufnahme Abb. 5 
zustande. Man sieht darauf deutlich das erste 


Abb. 5. 
Aufnahme einer Trépfchenbahn. 


2 Es wurde eine gewélbte Grenzflache gewahlt, weil bei 
einer ebenen die Flachenteile, die vor und hinter der 
Objektebene liegen, bei dem intensiven Gegenlicht das Bild 
teilweise tiberstrahlen wiirden. 
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Abb. 6. Berechnete Schnel- Abb. 7. Berechnetes Pha- 
leprofile (Kurven) senprofil (Kurve) 
und _ eingetragene und  eingetragene 
MeBwerte fiir ver- MeBwerte fiir die 
schiedene Frequen- Frequenz 60 Hz; 


zen. © linke Schatten- 
grenze, 


x rechte Schatten- 
grenze. 


Maximum der Geschwindigkeitsamplitude, wie 
es nach Abb. 2 zu erwarten ist. 

Solche Bilder wurden fiir die Frequenzen 30, 
60, 100 und 200 Hz aufgenommen und ausge- 
messen. Abb. 6 zeigt die Ergebnisse. 

In Abb. 5 ist innerhalb der Trépfchenbahn ein 


dunkler Streifen zu sehen. Er dient der Bestati- 
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gung der G1.(5). Der oben erwahnte Synchron- 
motor wird aus demselben elektrischen Generator 
wie der Lautsprecher betrieben; die Frequenz, 
mit der der Lichtstrahl unterbrochen wird, 
stimmt daher mit der Schallfrequenz iiberein. Da 
der Facher in die Objektebene abgebildet wird, 
ist die Trépfchenbahn an einer Stelle scharf 
unterbrochen. Weil die der Grenzflache nahen 
Luftschichten gegen die entfernteren phasenver- 
schoben schwingen (vgl. GI.(5)), verschiebt sich 
der dunkle Streifen in der Nahe der Grenzflache 
etwas zur Seite. Die Auswertung einer Phasen- 
messung zeigt Abb. 7. 

Die beschriebene Methode der Sichtbharmachung 
und Ausmessung der Schallausschlage (bzw. 
Schallschnellen) eignet sich auch gut dazu, in 
engen Kapillaren, deren lichte Weite mit der 
Zahigkeitswellenlange vergleichbar ist, die Schall- 
dampfung zu untersuchen. Uber diese Versuche 
wird in einer spateren Arbeit berichtet werden. 


(Eingegangen am 27. August 1952.) 
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ON A VARIATIONAL PRINCIPLE IN ACOUSTICS* 


by OSMAN K. MAWARDI 


Acoustics Laboratory, Massachusetts Institute of Technology, Cambridge 39, Massachusetts, U.S. A. 


Summary 


The present paper describes an investigation of the propagation of sound inside an acoustic 
labyrinth. The latter consists of a tube of rectangular cross-section interrupted at periodic inter- 
vals by thin partitions with staggered openings. By virtue of the expected periodicity of the wave 
function in the labyrinth, FLOQUET’s theorem can be employed to determine the velocity potential 
anywhere in the structure as a function of the potential in the space between any two consecutive 
partitions. This property is used in developing a variational principle which yields the propaga- 
tion constant of acoustic waves travelling down the labyrinth. The phase velocity and attenuation 
for the simpler type of waves propagated have been evaluated by the previous principle. The 
accuracy of the method is believed to be quite satisfactory. 


Sommaire 


On donne dans le présent article les résultats d’une étude sur la propagation du son dans un 
labyrinthe. Celui-ci est un tube de section rectangulaire partagé a intervalles constants par des 
cloisons minces ayant des ouvertures en chicane. Comme la fonction d’onde doit étre périodique 
dans le labyrinthe, on peut utiliser un théoréme de FLOQUET pour déterminer le potentiel des vitesses 
en tout point de la structure en fonction du potentiel dans l’espace compris entre deux cloisons 
successives quelconques. On utilise cette propriété pour établir un principe de variation qui 
fournit la constante de propagation d’une onde acoustique parcourant le labyrinthe. On a égale- 
ment évalué, par application de ce principe, la vitesse de phase et l’affaiblissement dans le cas 
d’un type simple d’onde se propageant dans le systéme. On estime que la précision de cette 
méthode est tout a fait satisfaisante. 


* This work has been supported in part by the Office of Naval Research, U.S. A.—Presented 
at the Eighth International Congress on Theoretical and Applied Mechanics, Istanbul, Turkey, 1952. 
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Zusammenfassung 
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Die vorliegende Arbeit beschreibt eine Untersuchung iiber die Schallausbreitung in einem 
akustischen Labyrinth. Dieses besteht aus einem Rohr von rechteckigem Querschnitt, das in 
periodischen Abstiinden von diinnen Zwischenwinden mit gegeneinander versetzten Offnungen 
unterteilt wird. Mit Hilfe der erwarteten Periodizitat der Wellenfunktion im Labyrinth kann das 
Theorem von FLOQUET benutzt werden, um das Geschwindigkeitspotential an einer beliebigen 
Stelle des Rohres aus dem Potential im Raum zwichen zwei aufeinander folgenden Trennwanden 
zu bestimmen. Diese Eigenschaft wird benutzt, um ein Variationsprinzip zu entwickeln, mit dem 
die Fortpflanzungskonstante der akustischen Wellen, die durch das Labyrinth wandern, sich 
bestimmen laBt. Die Phasengeschwindigkeit und die Dampfung fiir die einfacheren Wellentypen 
wurden nach diesem Prinzip berechnet. Die Genauigkeit der Methode scheint befriedigend zu sein. 


1. Introduction 


In several acoustic installations, constrictions 
are introduced inside ducts in order to increase 
the attenuation of sound waves travelling down 
the duct. A survey [1] of the technical literature 
reveals that the effect of discontinuous constric- 
tions on the propagation of sound waves has been 
extensively investigated. The mathematical 
methods used in these studies have borrowed 
heavily from analogous problems in electro- 
magnetic waveguides [2]. An exact solution for 
the propagation through discontinuous constric- 
tions, however, is difficult and often intractable. 
As a result, one has to resort to approximate 
methods. Of these, the variational techniques ap- 
pear to be the most fruitful, judging from the 
wide variety of cases which have been investi- 
gated by their means. 

The presence of periodic discontinuities in the 
duct adds appreciably to the complexity of the 
situation. Information on this class of problem is 
scanty and only a few isolated instances have 
been treated [3]. A renewed interest in the under- 
standing of the propagation down structures with 
periodic discontinuities was aroused, following 
the encouraging preliminary results of the use of 
an acoustic labyrinth [4] (Fig. 1) as an acoustic 
termination. 

The purpose of the present paper is to work out 
a generalized theory for the propagation of waves 
down the labyrinthlike structure of Fig. 1. The 
results are presented in the form of a propagation 
constant derived from a variational principle 
developed herein. 


2. Formulation of the problem 


A well known proposition in acoustics states 
that the velocity potential inside a region can be 
determined as a function of the value of the poten- 
tial or its gradient at the boundary. This func- 
tional relation follows from a direct application of 
Green’s second identity which leads to: 


P(t) =| (E(t) OnP(r') P(r’) enG(e,r’)]dS (1) 


where G(r,r’) is the Green’s function for the 


domain under consideration, @(r) is the velocity 
potential and r, r’ are position vectors referred to 
some origin 0, Choosing a Green’s function which 
has a vanishing gradient at the boundary, the 
above expression becomes: 


P(t) —[[G(r,r’) nD(e’)] AS. (2) 


S 


, 1el2 @ | 
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Fig. 1. Schematic diagram for labyrinth structure. 


Let a coordinate system x, y, 2 be defined (Fig. 1) 
such that the z-axis coincides with the “axis” of 
the labyrinth. Let it be assumed that the prevail- 
ing symmetry of the system will make the motion 
independent of the y-coordinate. For any cell n, 
enclosed between the two parallel planes z=nl 
and z=(n+1)l, it is found by means of eq.2: 


o (x, z) = tes (x, 2; bayer) np” (x, z') |e =n dS + 


Ss 


+ [[G (x, 23x',2") On O” (x’, 2’)]y— nant dS, (3) 


Ss 
The superscript n refers to quantities related 
to the domain nl <z, 2’<(n+1)l; 0Sx, x’ Say. 
Only one half of the labyrinth need be considered, 


by virtue of its axial symmetry. A similar relation 
derived for the cell (n+ 1) is, 


OO 8) (4) 


oe ‘i fers (x, 25 Si 4% ) eno" i (a, 2')|y =(n+l1)jl dS + 
S 


a f [G"*” (x, 25x',2') eno" +1) (x’,2')] = o so dS, 
5 


The domain of the cell (n+ 1) extends over 
(n+ l)IlSz, 2’<(n+2)l; OSx, x’Sa,. 


The labyrinth walls being rigid, the normal velo- 
city will vanish at the boundaries. By virtue of its 
symmetry, the x-component of the velocity will 
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vanish on the axis. The boundary conditions will 
further require the velocity to be continuous 
across the openings, i.e., for any n, 


8.0" [x, (n+1) 1] = 0.0" [x,(n4+1) 1] ay 
or 6.0") (x,nl) = 0.0" —" (x, nl) 0S «Sa, 


according to the case on hand. 

The propagation problem will be formally solv- 
ed, once the potential function is determined from 
the set of eqs. 3 and 4 in a manner compatible 
' with eq. 5. 


=%* Sa, 


(5) 


3. Construction of Green’s function 


The Green’s function for the cell n satisfies the 
differential equation 


2 
iS +t) 


6” (x, 23x’, 2’) =—3(x—x’) 8(2—2’) 


Ox? az” 
(6) 


subject to the conditions, 


OnG=0 at x=0, a, and at z=nl, z=(n+ 1). 
Let 6” = (7) 

es ole: pales, os ( = cos _ G" (z,2’), 
Whereoy7— io — ms Orn —0, 0 ==. 


A substitution of eq. 7 in eq. 6 shows that 
Ci mar 
24” (n) UNMEI OE Rye 
pares ) | es (z,z') = —S8(z—z’). (8) 


3 
The general solution of the above equation is (9) 


Go 


2 
where Ke =f — (re) and nlXz,z2’<(n+1)l. 
3 
Using the boundary conditions, the ratio of the 
arbitrary constants can be determined and it is 
found 
A sin Ky nl 


B cosKmnl’ (20) 


From an integration of both sides of eq. 8, it 
follows that 


z= 7 — 


0.6m (%2') 7) =—1 (11) 
cos (Kn nl) 
whence B= Re KT (12) 


The general expression for the Green’s function 


is finally found to be 
os (™™* 
a3 


ae yee 
a3 Kmsin Km if 


as 
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(Asin Km2'+ Boos Knz') cos Km [(n+1)1—2] \ 2° <2 
cos K,,[(n +1)1 —z'] (Asin Knz+ Boos Knz) Ja<z' 


oy a ada m(nl—2')] 2'<z 
‘eos | Kn [(n +1) 1—2z']} cos [Km (nl — 2)] eae 
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4. Derivation of variational principle 


The mathematical procedure in solving the 


' problem on hand is considerably simplified by 


introducing the concept of the propagation con- 
stant. Since the higher modes are expected to be 
rapidly attenuated it will be assumed that, 


D(x,nl) = ®[x,(n+2)l]er™", 


4P ifig ae 1 == Oy Woes (is) 


The previous expression (eq. 14) can be considered 
as following from a direct application of FLo- 
QUET’s theorem applied to the periodic structure 
studied. The constant « is in general a complex 
quantity, the real part of which vanishes when 
no attenuation is present. The real part cannot be 
negative when no sources are included in the 
labyrinth. 
The two eqs. 3 and 4 can be rewritten as, 


op” (x, 2) 
— [6 (x,252',nl) a, B(x’, 2’) 


0 


dx’ + 


2 — nl 


+[o" [x, 25x’, (n+ 1) T] Oy eo” (x’,2') dx’, 

2 Z#=(n+1)l 
(15) 

D+ (x,2) = fears, (n+ 1)1]- 
oy OO (a 2) dx’ = 
z= (n+ 1)1 
+ [6 fees’, (0 +2)1)- 
0 
a, pl" ote V(x’, 2’) dx’. 


Y= (n+21 

The negative sign appearing in front of some 
of the integrals follows from the accepted con- 
vention of letting the positive direction of the 
gradient be along the outwards normal. Since the 
velocity is continuous across the opening, then 
matching the velocity deduced from the two 
expressions (eq. 15) yields, 


20,0" (S532) ve ee 


a (16) 
= (1-6) f 2.6 (x,252',2') 0 O (x'2))) de’, 
o Sugiets jl 
a = 0S... 
(13) 


nlsSz, 2’<(n+1)l. 
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The last expression made use of the properties: 


Oz eo (4.25 Ste z')|s =(n+1)l — 
2 —nl 


=o) (ee Y(x,23x',2’)],— ee 
Pe esayt 


a if 
Oz ee (Ke 25 Ay | Pe +1)l=27 = ( ) 


4 n+ 1) 
= 0, [6° * jz; x".2) nee - 
By a similar line of argument a similar relation 


for the velocity distribution across the opening 
at z=nil, yields 


20: po”) (a.2)e a i 


= (1— er) [aC x,252',2' 6, @"")(x’, 2’) 


2=(n+1)l 
a z=nl 


VS4=c7 


For simplicity, let the following notation be intro- 
duced 

a.B" (x,2)s-n1=fi” (x), 

a.B” (x,2)s-401 =f2 (%)s 


then the question to be considered deals with the 
solution of the simultaneous integral equations, 


afi (x) = 


(19) 


(20) 
=(-l+e ) fee” (x,25%',2') |v sna” (x) dx’. 


Elimination of f\"(x) from these two relations 


yields, 
f= —F (Loo) 
4 
ff 16 (x,252', z'\]y=nl x 
et es ame) 
x Oz iC (xo phe aes pg) =(n+1)1° 
- fs” (x*) dx’ dx*(1—e°*'). 
Denoting 
fe (oe (Geer ea 2’) |y =nl 
S =(n+1)l 
‘ a:[C™ (x', 2; x*,2*)].*— aR 11 dx’ 


FY | 
by K"(x,x*) which by virtue of its construction 
is a positive definite function [5], then the above 
expression (eq. 21) can be rewritten, 
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fr” ( 


, (22) 
(cosh 2x1 — 1) [K" (x, x*) fs" (x*) dx*. 


a, 


x) 
1 
ye 
Inspection of the above relation (eq. 22) reveals 
that the propagation constant can be estimated 
by a variational formula. The eq. 22 is a canon- 
ical FREDHOLM integral equation of the second 
kind, and the unknown constant plays the part 
of the proper value of the equation. Thus multi- 
plying both sides of eq. 22 by f{"(x) and inte- 


grating between the limits a,—a,, we find 


1 
1 
9 (cosh 2x1 — 1) 
Jf 2° (2) K" (a,x) fo (x4) dcdx* 


The above expression is identified with the ratio 
of two definite positive quantities and as a result 
(cosh 2 x1 — 1) is stationary [6] for first order varia- 
tion of the function f;" (x). 


5. Evaluation of stationary quantity 


In order to determine the quantities to the left 
hand side of eq. 23, suitable trial functions f{"(x) 
have to be evaluated. The static solution is a 
reasonable guess. The latter is analogous to the 
determination of the lines of current in a conduct- 
ing plate ABCD (Fig. 2), the electrodes of which 
are at PB and DQ. This problem has been solved 
exactly by Moutron [7], who gave the distri- 
bution in terms of elliptic functions. A close ap- 


— 


D 


B Cc 


Fig. 2. Idealized conducting strip used to calculate the dis- 
tribution of velocity in the aperture. 


proximation found by assuming the constrictions 
far apart (i.e, AD>DC) yields results more 
amenable to easy computations. This problem 
also had been previously worked [8]. The complex 
potential in the vicinity of the aperture say at 
z=(n+1)l(a,<x<a,) is given by the results of 
ref. [8] as 
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sin ee 
w = —  arcsin ee a (24) 
a. a , 
a sin $93 __°2) 
Le 


where Z is the complex number (x +jy). When the 
arbitrary constant A is made equal to a, then the 
flux of fluid through the aperture is given by 


as 


rf fk (22) dx = ag. (25) 
The kernel is found to be: 
Kix. si Bi: cos? kl + 
a3 (26) 


* 
a . 2a x eX 

+ 3 cos K,l sin — (cos — + cos ——}) cosk 1+ 
4 as as as 


I 2 s 
+= cos? K,l(a,+ #8 sin "| ( cos cos ) +... 
2 27 as, as a 


Substitution of eg.24 and eq. 25 into eq. 23 
yields for the propagation constant 
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Attenuation in Nepers per section 


Kl 


Fig. 4. Real part of propagation constant corresponding to 
Fig. 3. Only one attenuating band has been com- 
puted. 


quite rapidly, only the first few terms have been 
retained. 

In order to investigate the attenuation bands 
for the labyrinth C/2D, a typical case has 
been drawn in Fig. 3. Attenuation is present 
whenever 0>C/2D>2. The corresponding pro- 


cosh 2a1 = 
3/2 (1%?) + 3/2xsin™™ 
a3 a3 
Lia a Ma,/a, . ma It joe al Qna,\ . ,na ars ror 
2 |— coskl — sin? —?! 3 sin cos K,lcoskl4 sin chee rae eicos? Kala. ; 
a; Xt 3 2\a, 22 a, 2a; 


which can be written formally, 


cosh 201 = 1—C/2D. 


<a eb) 


2 4 6 


kl 


Fig. 3. C/2D as a function of kl for the case of a,/a,; = 0.8, 
a,/a,;=0.5 and l/a,; = 3/z. Attenuation occurs when- 
ever 0 >C/2D >2. 


The propagation constant « then becomes a pure 
imaginary whenever | 1 — C/2 D| <1; it, however, 
will have a real positive part indicating attenu- 
ation in the propagated wave, whenever 
|1- C/2D|> +1. Since the series of D converges 


pagation constant has been computed and is 
shown in Fig. 4. 


(Received 6th November, 1952.) 
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ACOUSTIC IMPEDANCE MEASUREMENT 
BY THE TRANSMISSION-CHARACTERISTIC METHOD 
by A. F. B. NICKSON and R.W. MUNCEY 


Division of Building Research, Commonwealth Scientific and Industrial Research Organization, Australia 


Summary 


The transmission-characteristic method for measuring acoustic impedances has been developed 
to the stage where it is an accurate, simple and versatile method. The most suitable chamber 
dimensions are derived, construction of chambers is described and the impedance of a selected 
material at various angles of incidence is given for a wide range of frequencies. 


Sommaire 


On a perfectionné la méthode de mesure de l’impédance acoustique, dite méthode de la carac- 
téristique de transmission, de fagon a la rendre précise, simple et universelle. 

On a calculé les dimensions les plus convenables 4 donner a la chambre, on décrit le mode de 
réalisation de cette chambre et on donne, pour une gamme étendue de fréquences, l’impédance 
d’un matériau choisi, mesurée sous des angles d’incidence différents. 


Zusammenfassung 


Die Methode der ,, Ubertragungs-Charakteristik* zur Messung von akustischen Impedanzen ist 
zu einer genauen, einfachen und vielseitigen Methode entwickelt worden. Die giinstigsten Kammer- 
abmessungen werden abgeleitet, der Bau der Kammern wird beschrieben und die Impedanz eines 
ausgewahlten Materials wird fiir verschiedene Einfallswinkel und fiir einen weiten Frequenz- 


bereich mitgeteilt. 


1. Introduction 


The classical methods for measuring the acous- 
tic properties of bounding surfaces employ the 
reverberation chamber which yields the absorp- 
tion coefficient and the impedance tube which 
measures the specific normal acoustic impedance. 
Both methods have serious shortcomings. In a 
reverberation chamber it is difficult to achieve 
satisfactory diffusion for frequencies below 500 
to 1000 c/s, the absorption coefficient is some- 
what dependent on the size and shape of the 
chamber in which the test is made and in actual 
use the coefficients seem to be too high. The short- 
comings of the impedance tube are that it oper- 
ates at normal incidence only and has artificial 
mounting conditions. 

Consequently it is strange that the transmis- 
sion-characteristic (t.c.) method first proposed by 
WeENTE [1] and used by Harris [2] and others 
has not been widely adopted; with it measure- 
ments can be made at definite angles of incidence, 
mounting conditions can be close to standard and 
the use of low frequencies is comparatively simple. 
High accuracy is required only in measuring audio 
frequencies. As this can be done by comparison 
methods, it is not difficult. The Division of Build- 
ing Research has been investigating the t.c. 
method and is finding it most useful. 


Harris has shown that the impedance of a 
material can be calculated from the change of 
bandwidth and resonant frequency of the normal 
modes of vibration of a rectangular chamber 
when the material is substituted for one wall of 
the chamber. His equations are given below. His 
analysis is not rigorous in that he does not actu- 
ally consider the case of forced vibrations; never- 
theless the equations he derives are the same as 
those obtained by a more rigorous treatment. 
Consider a (bare) rectangular chamber whose 
dimensions are l,», ly, and l., and suppose that 
material is introduced to cover the y-z wall so 
that the new x-dimension is ],. Then, use the 
following symbols: 
vy = the frequency of the mode of vibration with 
no material in the chamber, 

ve = the frequency of the mode of vibration with 
the material in the chamber, 

ky = the bandwidth with no material in the 
chamber (taken as half the frequency dif- 
ference between frequencies at which the 
response has dropped 3dB from maximum), 

k. = the bandwidth with the material in the 
chamber, 

o = the density of the air, 

c¢ = the velocity of sound in air at the particular 
conditions, 
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Nx, Ny, Mz = the number of nodal planes in the x, 
y and 2z directions, 


Ox, Gy, 6: = coefficients of value % or unity de-. 
pending on whether the correspond- 
ing Nx, Ny, Nz are zero or non-zero, 

A = the wave length of the sound, 

Weal Vey hy =... CbC., 

Coc = the impedance of the material, 

Z =ocyexpj®=the impedance of the 

chamber wall. 
Harris shows that the impedance of the mate- 
rial measured is given by 
. 2 coth(—jzyx) 
aT iy 


where the real and imaginary parts of y? are 


) 7 Z 2 reas 2 2G5 Ey sin @ 
Re (2) =n? (I) Ah (mom). See a 


9 


ES G Tt 4 
4vl, 
Im (7?) = ey 
axl bles ky 
: ke k A TF eM Ipr: 
( oy 2 (ox Lye Ib + Oy Lp Les + Oz lxp ‘| 


2. The chambers 


The transmission-characteristic method uses 
normal modes of vibration of the chamber which 
are separated in frequency (or which can be sepa- 
rated by microphone placement). It is desirable 
therefore to proportion the room so that the 
separation in frequency of the lower modes is as 
uniform as possible. 

The frequencies of the normal modes of vibra- 
tion are given by 


Ley V ayy ley Lin\* 2 
2 G =Vni felis Ny + (7) Nz 


By writing (lx»/ly)? =D and (Ixs/ls)?=E it can 
readily be seen that the two modes of vibration, 
(P,, Q;, R,) and (P,, Qo, R.), will have the same 
frequency if and only if 
Pi + DQi + ERi = P; + DQ: + ER}. 

The locus of points satisfying this condition is 
a straight line on the D, E plane. Fig. 1 shows 
all such loci for which the values of P,, Q,, ete. 
are less than three. Now it is desirable that the 
dimensional ratio chosen for the chambers should 
be such that the corresponding (D, E) points do 
not fall upon the loci and, in fact, the distance 
away from a line is some measure of the separa- 
tion of the two modes in question. Hence it is 
necessary to fix the dimensions so that the cor- 
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responding point falls in one of the spaces of the 
diagram. 

The examination can be made as follows. The 
ratio 1:D:E is chosen to fall in an open space of 
the diagram (e.g. H, J, K) and the modal fre- 
quencies calculated based upon the lowest modal 
frequency equal to unity. The frequency ratios 
corresponding to three areas H, J and K are given 
in Table I. These and other series were examined 
for uniformity of spacing and separation and the 
ratio 1:0.785:0.5 chosen as the most satisfactory. 


— Value of E 


100 
100 
ASS 

0.75 =k { 
es W. / 
= SD & 
3 KKIKIG i 
S 050 AZ 
| Hycboe YY 


Fig. 1. 
Loci on the D-E plane on 
which two or more normal 
modes (N=0, 1 or 2) have 
the same frequency. H, J, 
K = areas in which desirable 
dimension ratios may lie. 


Table I 


Frequencies of normal modes for various chamber 
dimensional ratios 


1:0.96:0.43 | 1:0.785:0.5 | 1:0.52:0.43 


Dimension 
ratio 


Modal 1.0.0 1.00 1.00 

frequencies} 0.1.0 1.25 1.10 
0.0.1 1.4] 52 
1.1.0 1.60 1.49 
1.0.1 1.73 1.82 
0.1.1 1.89 1.88 
2.0.0 2.00 2.00 
Tava 2.13 2.13 
2.1.0 2.36 2.28 
2.0.1 2.45 2.51 
0.2.0 2.50 2.20 
1.2.0 2.69 2.42 
Pee | 2.74 
0.0.2 2.82 
0.2.1 2.87 


Having fixed the dimension ratio it is necessary 
to decide the largest dimension of each chamber. 
An overlap of frequency range is desirable and a 
ratio of about two to one appears reasonable. The 
useful frequency range can be covered by six 
sizes as follows: 


194, 


Lowest mode 


frequency 100 200 400 800 1600 3200 e/s 


Greatest chamber 
dimension 


1.73 0.864 0.432 0.216 0.108 0.054 m 


Chambers have been constructed so that at 
least two are of the same material and have the 
same wall thickness making them thermally simi- 
lar. When such a pair is placed in the same ther- 
mal environment the internal temperature (and 
consequently the wave velocity) of the two will 
drift equally. This drift is reduced in the labora- 
tory as the space in which the chambers are 
placed is insulated to lessen the rate of tempera- 
ture change. 

The two largest chambers have been construct- 
ed in 10.2 cm concrete blocks hard-plaster- 
rendered inside, those with greatest dimensions 
0.864 m, 0.432 m and 0.216 m have been made 
from 2.54 cm (nominal) wood lined with 0.16 cm 
aluminium sheet, and the three smallest are con- 
structed in 1.48 cm brass. The masonry chambers 
have concrete doors with refrigerator type seals 
and locks and the samples are placed inside; the re- 
mainder have one open side which can be cover- 
ed by the lid or by the sample and in this way 
measurements can be made on in situ materials. 

In the masonry and wooden boxes, sound is 
produced by a small loudspeaker mounted in one 
corner, and in the diagonally opposite corner a 
small microphone connected to a vacuum tube 
voltmeter measures the amplitude of the sound. 
In the brass boxes sound is introduced from out- 
side by a tube stuffed with wires and is detected 
by a microphone also outside and connected to 
the box with a similar tube. As indicative of the 
wall hardness it is sufficient to say that in the box 
1.73 x 1.37 x 1.22 m? (volume of 2.89 m3) the 
R.T. is 13.5 s at 100 c/s and in the brass box 
10.8 x 8.58 x 7.63 em? (volume of 706 cm?) the 
R.T. is 0.5 s at 1600 c/s. These decays correspond 
to Sabine absorption coefficients of 0.004 and 
0.007 respectively. Consequently in the present 
measurements the walls have been taken as being 
infinitely hard and no correction is applied in the 
second equation. For experimental purposes the 
expressions for Re(z”) and Im(y?) of Section 1 
may be approximated to 


2 2 ‘LE ‘ Qn, Av 
meets (i) ery, (1) 
Hallie 
and = Im (7?) = a (ke — mks) (2) 
where Ay =e. 
and mt o Ox ly La 


2 (oxlyoles + dy Lgles + ole lys) ” 
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3. Experimental procedure 


The equations (1) and (2) involve accurate 
determination of small frequency increments. The 
required accuracy is obtained by measuring each 
frequency in terms of a standard frequency of 
good stability. A stable oscillator, the frequency 
of which can be varied, is used to locate the par- 
ticular frequency to be measured. A second oscil- 
lator is adjusted to a frequency close to that of 
the first oscillator (the second frequency being 
known exactly in terms of the standard crystal 
frequency) and the frequency difference between 
the two oscillators measured by means of the 
Lissajous figure formed on the screen of a cathode 
ray tube when the two frequencies are applied 
to the deflecting plates. In this manner the fre- 
quency of the variable oscillator is determined to 
a high degree of accuracy, the estimated error 
being less than one part in 10 000. The equipment 
is described in detail elsewhere [3]. 


4. Calculations 


The calculations required can be shortened by 
using conformal transformations in the manner 
shown by Morse and Bo tr [4] who give several 
maps covering a range of impedances for the 
transformation of y? to In(¢/n). For many ex- 
perimental measurements the maps are not suf- 
ficiently contoured or extensive and failing con- 
siderable calculation and re-drawing the follow- 
ing procedure was adopted: 

sf coth (— jz x) coth W 
f=. —— / = — aN: ay 


where W ex jy. 


coth W/W can be expanded as terms of a series 
containing one set of real terms and one set of 
imaginary terms. When this is done it is found 
that the first two terms in each set are sufficient 
to cover most of the acoustic impedances in which 
interest is shown. Equating the real and imagi- 
nary parts of ¢ with the real and imaginary parts 
of the expansion multiplied by the quantity 7, 
gives two equations for the two parts of ¢ in terms 
of the real and imaginary parts of 7. (These latter 
are determined from experiment.) The results 
obtained in this manner can be expressed as 
follows: 


(I).n==.0: 
R _ Cc ad (ke — mk») 
oc =A, (ke — mk»)? + (v5 — ve)?” 
i a ee (seek 
oc 3 4a, (ke —mki)?+ (v— v%)?’ 


ACUSTICA 
Vol. 3 (1953) 


(II) nx=1,2...1.= 1. (box with removable side) : 


Ty pega uneaaae CRAEL 
oc §= Axl, (ke — mks)? + (vs — »)?’ 
xX c — (v¥ — ve) 


oc Qal, (ke — mks)? + (v% — ve)?’ 
(III) n, = 1,2... 
side) : 
R fa (s f (ke = mks) 
oc al, (ke — mk»)? + [(v — ve + «) (Le/lev)?]?’ 
4 ce — (% — Ve + x) (Lx /Lee)? : 
oc 2al, (ke — mks)? + [(v% — ve + x) (Ix /leb)?]2” 


1. + 1.» (box without removable 


2 
nv 


wa (1 —le/l). 


where i 


The accuracy of these relations is such that the 
maximum errors are less than 5 per cent for 
Waies ue Vb— Ve ue ft ote ey ke —mkp O04 
v v Zs 
x Nx 
with 1, /Lep == 0.95. 

Note in (I) and (II) that a material of half the 
R and X impedance is necessary to give the same 
bandwidth and frequency shift for a grazing 
wave as for a normal or oblique wave, which 
agrees with the values used for o (o=% for 
n=0, o=1 for n=1) since usually 27/3 « X/oc. 


5. Bare wall impedances 


The value of the impedance of the test specimen 
has been expressed in terms of the frequencies of 
the normal modes of vibration and the resonance 
bandwidths associated with those modes by 
equations (1) and (2). Each equation contains a 
correction term for the impedance of the wall op- 
posite the sample, the effect of the other bare 
walls being eliminated by the analysis. In the 
first of these equations the correction term can 
easily be computed but in the second a knowledge 
of the bare wall impedance is necessary. 

Harris [2] points out that the impedance 
could be calculated from a knowledge of the 
sound velocity and the room dimensions, but re- 
jects this method because of the large error likely 
to be introduced in the impedance by small errors 
in the sound velocity. He then proposes that the 
pressure distribution in the chamber be explored 
and, assuming equal wall impedances for the case 
when nx>1, shows how the impedance may be 
evaluated. This technique cannot however be 
applied to the case where ny=1 because the 
minimum pressure is theoretically zero and con- 
sequently the impedance of any wall of a chamber 
cannot be determined below twice the lowest 
modal frequency. Further, the wall impedances 
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have been found to be so high that it is doubtful 
if measured minima would give as good an accura- 
cy as the calculations using the velocity of sound. 

It was decided therefore to adopt the method 
using the velocity of sound and to assume its 
value to be 340.5 m/s at 15°C with a temperature 
coefficient of 0.602 m/s per degree centigrade. 
These figures were based on evidence from various 
sources, notably ABBEY and Bartow [5] and 
Harpy, TELFAIR and PIELEMEIER [6]. This value 
may be expected to be within +0.2 m/s and 
such an error has no significant effect upon the 
order of the calculated impedance. 

Morse and Bo tt [4] have shown that the fre- 
quency factor for a standing wave when all the 
walls are fairly hard has the form exp (—j@nt —knt) 
where wn and k, can be written 


2 2 RE Spar ari) 2 
On =T7C ve s a Py = a ee ane 
- i P 
c i Mxkx 400 My Ky An [lz Kz 
kn ee a — 
4 Vx L Ny EF Vez Ll, 


and the various y’s and x’s are the respective 
wave number and attenuation parameters. 
A series expansion for w, and k, gives the value 


2 2 20 
[x —k, =n? — ae (s + Sp) 


ATC Mx Kx 


and Bere Ao (g, + 8) 
priOh anGh 
where s=7 sin Orgs 7 08 ®. 


The approximate value of k, for the damping 
constant when all walls are fairly hard is therefore 


kn —- 


¢ [40x (Bat Bex) , 40y(Byi + By2) | 402 (8a + 822) 
4 I. i l. 


or if values of g are the same for each surface of 
the box, k, is found to be equal to 


cos D Ox Oy Os 
Re Laie where Bearcat Fok Lr: 
In a similar manner a value may be obtained for 
Ons 
2 20x Nx 


Nx (Sx af Sx2) 
2 2.2 It 
On = 70" c* | — om Sn a 


Remembering that (2) eal ls ds Sp onl 

Cc, Eolebs 
and calling this frequency ») as being the value 
if the walls are infinitely hard, the value for w is 
found to be 
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c 7 C 


wo. =7 ce iI 88% Zs K| 


and q@- is the value of frequency measured with 
fairly hard walls. This equation may be put into 
the form 


os Ye a PKOAK i 2) 
wi le ee CAL 


From equations (1) and (2), the value of Z= 
ocy + expj @® can be obtained from the solutions 


es ROIS f ee (2 — a) Le 


Vol kn 


Vbi TU 


4 (ve — vin) 
V1 kn 


. 


@ = arc tan 


Measurements made on the softest of the three 
wall types used (the wooden surface lined with 
sheet metal) resulted in a specific impedance 
120 /—60 at 400 c/s. This gives a value of sin O/y 
as 0.0065 that is small compared with the value of 


ane” and may therefore be neglected. 
v 


6. Results of experiment 


To illustrate the use of the t.c. method for 
measuring the acoustic impedance of a material, 
it was decided to measure one material over a 
large frequency range, and for as many angles of 
incidence as could be conveniently measured in 
any one box. As the normal modes in any box 
come closer together and the bandwidths become 
wider as the frequency increases, there is a limit 
to the range of frequencies that can be measured 
in any box. For the material measured this ratio 
was about 3:1 in the larger boxes and 2:1 in the 
second smallest box. In the smallest box only a 
small range was covered because of inefficiency 
of the loudspeaker and microphone units when 
used with tubes filled with fine wires. The material 
for which the measurements were made was a 
soft fibreboard, 1.27 cm thick, weighing 0.38 g/em?. 
Samples of this material from a standard batch 
were cut and glued into position in the boxes with 
a cold water glue. The samples were allowed to 
remain for five days before measurements were 
made. The experimental procedure was to meas- 
ure the frequencies and bandwidths of modes in 
one box before the sample was glued in, referring 
at intervals to any one selected mode in a thermal- 
ly similar chamber. In this way all the modal 
frequencies of the chamber could be obtained as 
a ratio using the arbitrary selected mode which 
was kept under observation. Five days later when 
the glued sample was ready to be measured, the 
modal frequencies in the chamber were again 
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measured and the new bandwidths determined, 
referring occasionally to the initially selected 
mode in the other thermally similar chamber. As 
this would have changed from the previous oc- 
casion when measured, this value in conjunction 
with the ratios determined at the start will give 
the values of the bare box modes as they would 
be if there was no sample present. The new modal 
frequencies being then measured enable the fre- 
quency shifts to be caleulated. The bandwidths 
of these peaks are also measured. 

In some cases the peaks interfere with one an- 
other, and this difficulty can be overcome oc- 
casionally by using the decay of sound at the 
particular frequency being considered. The simple 
relation k,=6.91/T relates bandwidth to time of 
decay, and thus the bandwidth can be measured 
when it is not possible to reduce the peak inten- 
sity sufficiently on either side to enable k, to be 
measured in the normal way. A Briel and Kjer 
level recorder (type 2301) was found useful for this 
purpose although its resolving time was rather 
too long for some cases. Varying the position of 
the microphone also helps to operate modes when 
close together. 

Table II gives the values of R/oc and X/oc 
measured for the sample at various frequencies 
and angles of incidence for a number of different 
boxes. The greatest practical difficulty was that 
although the samples were cut from two sheets of 
a known batch, each sample had necessarily a 
different glued joint. All that could be done here 
was to mix and apply the glue in the same way 
each time, place weights on the surface and then 
measure the samples after five days when the 
joints could be considered reasonably similar. 

Table III gives the mean values of R/oc and 
X/oc when divided first into frequency groups 
100—200, 200—400 c/s, ete., and angular groups ~ 
0—30, 30—60, 60—90°. This enables an overall 
picture of the behaviour of the fibreboard for 
three different angles of incidence and a large 
range of frequencies to be obtained. These values 
are plotted in Fig. 2. 


7. Discussion of results 


To discuss the results profitably it is necessary 
to refer to the average results of the measure- 
ments shown in Table III and Fig. 2 rather than 
to the individual measurements in Table II. The 
important feature is that the acoustic impedance 
of the material depends on the angle of incidence 
of the sound wave. This is clearly shown in the 
figure, where it can be seen that the impedance 
for incident angles 0—30° is greater than for 
angles 30—60°, and the latter group has a greater 
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Mode Frequency R/oc X/oc 6 
Box 5 

001 1551 10.2 —11.5 81 
010 1951 11.2 — 8.59 83 
100 2210 6.71 — 8.42 0 
O11 2498 8.18 — 4.66 84 
101 2709 9.50 — 1.08 42 
110 2941 10.3 — 3.54 47 
002 3130 8.11 — 2.08 87 
OT: 3313 8.49 — 5.02 54 
012 3681 7.93 — 6.65 84 
102 3803 6.31 — 2.53 59 
020 3879 3.14 — 4,95 80 
021 4059 0.90 — 3.42 73 
112 4231 2.76 — 7.46 64 
Box 6—sample A 

001 3165 87 
010 3920 88 
100 4382 0 


Oll 4847 


Table III 


Frequency Mean Mean Range 
range [c/s] R/oc X/oc of angle 0 


Zz 


100—200 
200—400 
400—800 
800— 1600 
1600 —3200 
3200 — 64.00 


100—200 
200 —400 
400 —800 
800— 1600 
1600 — 3200 
3200—6400 


100—200 
200—400 
400—800 
800— 1600 
1600 —3200 
3200 —6400 
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Fig. 2. Frequency characteristic of 1.25 cm fibreboard at 
different angles of incidence. 
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impedance than for angles 60—90° over most of 
the frequency range from 100—5000 c/s. This 
means that for this sample the absorption co- 
BESS | will increase 
Re x2 

as the angle of incidence departs from the normal. 
This is to a certain extent a confirmation of the 
suggestion that the absorption coefficient would 
alter with angle of incidence of a material of 
normal impedance Z/oc in the manner of 
Zcos6— ocl? 
Zcos0+ oc|” 
mum at = arccosoc/Z and is zero when = 77/2. 
Although there are several modes in the boxes 
for which geometrically the angle of incidence is 
z/2, it is found in practice that the presence of 
the material modifies the wave behaviour in 
such a manner that the angle of incidence is not 
greater (in the present cases) than 88°, and would 
be less than this for a material of a smaller normal 
impedance. The manner in which the curves cross 
around 200 c/s when the impedance of the sample 
for angles 30—60° becomes smaller than for 
angles 60— 90° cannot be explained in this way but 
the effect is quite real for the material used. In 
the case of the lower frequencies, the above for- 
mula would not be satisfactory as a representa- 
tion of the variation of « with §. This method 
however gives the value of R/oc and X/oc for an 


efficient calculated from « = 


This expression has a maxi- 


There are several difficulties in using the t.c. 
method for measurement of acoustic impedance 
and some attention has to be paid to them. The 
obvious one is that only materials of medium 
impedance can be measured satisfactorily, be- 
cause after the first two or three modes the band- 
widths become so wide as to mask each other and 
make it difficult to determine the presence of peak 
frequencies. If this small number of modes is 
satisfactory, the range of incident angles is limit- 
ed to 0° or 90° and it is difficult to measure any 
other angles. The fibreboard in the last range has 
R/oc=5.8 and X/oc= —5.0 giving «=0.79 and 
this would seem to be about the limit of the meas- 
urement. It is also found that some of the reson- 
ance peaks are not symmetrical and this leads to 
difficulty. In cases where the decay could be 
measured with the sound-level recorder, the band- 
width so determined was generally found to be 
closer to the narrower half band than the wider. 
This was certainly so when one band was being 
interfered with by the next one. Moving the 
microphone position was sometimes found to be 
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helpful as this is equivalent in effect to suppress- 
ing the interfering mode. 

The third difficulty is due to the change in the 
velocity of sound with temperature. This is im- 
portant as small frequency shifts are often to be 
measured, and it has been overcome by using two 
thermally similar chambers in each measurement. 
It is important to watch that the frequency shifts 
measured are due solely to the sample and not 
partly to a change of temperature while the meas- 
urements are being made. 


8. Conclusions 


This investigation has shown that the trans- 
mission-characteristic method is a powerful tool 
for measuring the acoustic impedance of materials 
at various frequencies and various angles of inci- 
dence. Certain limitations are to be expected when 
measuring the softer materials and the method 
then requires modification. The obvious difficulty 
is due to increasing bandwidths while the modal 
frequencies are coming closer together, but this 
may be partly overcome by varying the position 
of the microphone in the box. The boxes may be 
portable and can be easily placed against mate- 
rials already erected for which measurements are 
required. The method provides an excellent means 
of comparing the impedance of various materials 
at different frequencies as required by the ex- 
perimental work on acoustic models. The fact 
that allowance can be made for the bare wall 
impedance means that for most measurements 
boxes made of wood lined with metal sheet will 
be satisfactory provided proper attention is given 
to the sealing of the box. 
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The mercury thermometer 


An instrument for the measurement of high-intensity 
ultrasonic waves 


by O. Linpstr6M 


Division of Physical Chemistry, 
The Royal Institute of Technology, Stockholm, Sweden 


Zusammenfassung 

Ein gewohnliches Quecksilber-Thermometer wird als 
Leistungsmesser fiir Ultraschall angewandt. Diese ein- 
fache Methode ist iiberraschend genau, der mittlere 
Fehler betriagt etwa ein Prozent. 


In the quantitative study of the chemical and biological 
effects of ultrasonic waves the determination and control of 
the sound intensity is a question of fundamental importance. 

The instrument readings on the power control panel of 
the high frequency generator, e.g. of plate voltage and 
plate current of the rectifier, yield accurate expressions for 
the electric power input. Under constant conditions there 
is in general a linear relationship between the electric power 
input and the acoustic output of the transducer. The quo- 
tient between the acoustic output and the electric input 
may be called the power conversion ratio, p.c.r. for brevity. 
Unfortunately the p.c.r. is not a stable apparatus constant 
which may be evaluated once for all but a rather changeable 
magnitude influenced by many factors difficult to control. 
Because of this variancy it is necessary in accurate work to 
make frequent determinations of the p.c.r. in order to 
secure a defined acoustic output. 

It should here be pointed out that the p.c.r. and the 
underlying linear relation refer to the primary acoustic out- 
put from the transducer. In general the sound intensity is 
evaluated by means of any instrument placed at some 
distance from the transducer. When the liquid in the inter- 
space between the transducer and the instrument is cavitat- 
ing the p.c.r. as determined in this way will decrease be- 
cause of the reflexion of sound at the gas bubbles formed in 
the cavities. Cavitation has therefore to be avoided, e.g. 
by working at low acoustic intensities, using degassed li- 
quids or increasing the hydrostatic pressure in the liquid. 
A convenient method is to arrange a fast circulation of the 
liquid so that the gas bubbles are swept out of the sound 
field before they have attained a disturbing magnitude. 
This was done in the experiments reported here. (At high 
rates of liquid circulation the transmitted power will be 
reduced due to an effect of the streaming liquid, therefore 
it is necessary to maintain constant conditions even with 
regard to the liquid circulation [1].) 

A great many methods have been designed for the deter- 
mination of the intensity of powerful ultrasonic waves in 
liquids. Most of them are based either on the sound pressure 
or on the heating action of the acoustic waves [2]. The for- 
mer methods, in which the sound pressure is measured by 
means of a suspended reflector, yield absolute intensity 
values. However, these methods are complicated because of 
the tendency of the reflector to remain fixed in certain pre- 
ferred positions, a phenomenon which is observable even in 
the case of a (slightly) conical reflector. The thermal 
methods do not easily give absolute values but are poten- 
tially more simple and rapid. An accurate thermal method, 
calibrated with the aid of an absolute sound pressure me- 
thod, therefore seems to be most suited for the continuous 
control of the p.c.r. during a series of experiments. It was 
then realized that an ordinary mercury thermometer might 
offer some advantages, the mercury bulb being sufficiently 


0. LINDSTROM: THE MERCURY THERMOMETER 


199 


large and the instrument itself a most reliable temperature 
indicator. The degree of reproducibility of such a thermo- 
meter method was evaluated in the following experiments. 

The most suitable position for the thermometer was 
found to be the one demonstrated in Fig. 1. It is easy to 
arrange the thermometer in this symmetrical position with 
the aid only of an ordinary laboratory stand, and the error 
in the position thus achieved is of the order of 1 mm. The 
important reason, however, for arranging the thermometer 
bulb in the acoustic fountain and not in the liquid bulk is 
that the surface of the fountain acts as a kind of curved 
reflector, which focusses the acoustic waves on the thermo- 
meter bulb so that a much larger temperature difference is 
obtained. 

In these experiments the 
liquid medium consisted of 
transformer oil, which was 
circulated and thermostated 
to 24.22 + 0.10°C. The con- 
stant oil temperature is a 
prerequisite in order to se- 
cure a fairly constant p.c.r. 
value, because of the tem- 
perature variation of the 
transducer’s resonance fre- 
quency. (Of course this does 
not imply that it is necessary 
to work with a thermostated 
liquid when applying this 
thermometer method.) The 
statistical variations in the 
acoustic energy (equal to 
+5% and due to fluctua- 
tions in the power supply) 
and the variation in the oil 
temperature made it necessary to evaluate the appropriate 
temperature difference as a mean value over a short period 
of time. Simultaneously the variations in the electric input 
were recorded by means of a printing potentiometer so that 
an accurate mean value for the electric input/could be com- 
puted. The small differences in the power input between the 
separate experiments were corrected for with the aid of the 
approximate linear relationship between the temperature 
difference and the power input. It has to be pointed out 
that the separate runs are independent of each other since 
the thermometer was taken away and the laboratory stand 
removed between each experiment. 

The results of some runs with four different thermo- 
meters are given in the Table. It is recognized that in spite of 
the rather simple procedure the standard deviation is suf- 
ficiently small, of the order of one per cent for the two more 
accurate thermometers A and D. This error could probably, 
if necessary, be reduced by means of some mechanical ar- 


Fig. 1. 
Position of the thermometer 
in the ultrasonic field. 
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Fig. 2. Temperature difference versus electric power input. 
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rangement permitting a more exact and defined thermome- 
ter setting. 

In Fig. 2 a typical curve for the temperature difference 
versus power input is reproduced. As the thermometer bulb 
at low intensities is partly situated in the air there is some 


Ther- | Seale | Hg bulb T ; Mean | Max. 
; 2 emperature diff. 

mo- div. d l value | error 

meter’ | [Clem] | _ "femal (°C) rc] | (%) 


9.35 9.12 
9.25 O31. 
7.94 TO 
7.96 8.15 
8.75 8.80 
8.66 8.43 
7.93 8.01 
8.20 8.06 


The teaching of acoustics 

A joint meeting of the Acoustics Group of the Physical 
Society and of the Education Group of the Institute of 
Physics was held at the office of the Institute in Belgrave 
Square, London, on the evening of the 19th November. The 
meeting, under the chairmanship of Mr. WEST, was very 
well attended and a live discussion followed the opening 
speeches, 

Dr. R. W. B. StepHEeNsS (Imperial College) in the open- 
ing speech stressed the study of vibrations as being fun- 
damental to physics as a whole and pointed out the close 
association of modern acoustical theory with that of other 
branches of physics, e.g. the similarity of DEBYE’s theory 
of specific heats with the theory of room acoustics. He sug- 
gested that the treatment of acoustics as a Cinderella sub- 
ject in schools and colleges arose from the crowded nature 
of the present syllabuses and the very small proportion of 
acoustic questions set in physics examination papers. The 
speaker then pointed out how acoustics permeated into the 
fields of aeronautics, light electrical engineering, rheology, 
medicine, etc., and examples quoted from this wide sphere 
of application could do much to make the teaching of the 
subject more alive and attractive. 

Mr. E. NIGHTINGALE, an experienced teacher and writer 
of school text-books, dealt with the teaching of acoustics in 
schools. He indicated, by means of demonstrations, how 
much of the earlier teaching can be illustrated experiment- 
ally with simple apparatus by the energetic and enthusiastic 
schoolmaster. A suggested list of topics for lecture demons- 
trations was presented by the speaker. 

Dr. E. G. RicHarpson (Durham University) dealt with 
the University aspect of acoustics instruction and research 
and pointed out that while in England there are no instit- 
utions or university departments solely devoted to the 
subject, yet advanced teaching exists in some five or six 
colleges and universities in the country. He pointed out the 
need for a research worker in acoustics to have undergone a 
balanced training in hydrodynamics and applied electronics. 

Mr. D. M. A. Mercer (Southampton University) gave 
tentative details of the acoustic syllabus for honours physics 
students at his University consisting of approximately 
7 lectures per year of a 3 year’s course. The speaker sug- 
gested that acoustics could be conveniently divided into 
3 sections, namely classical, modern and electro-acoustics, 
and he stressed the importance of electrical analogies and 
circuit theory in present-day acoustical practice. This point 
was also emphasised by other speakers, and also the need 
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uncertainty at the smallest values. This error vanishes 
when the acoustic fountain is high enough to cover the 
whole bulb. 

Longer continuous exposures of the order of thirty 
minutes or more must be avoided since the thermometer 
may then be damaged due to cavitation which 
may be set up in the mercury bulb. 

The writer wishes to thank Professor O. LAMM 
for his kind interest and encouragement. Thanks 
are also due to the Swedish State Council of 
Technical Research. 


Stockholm, 7th September, 1952. 
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for an adequate mathematical background for the ad- 
vanced worker in acoustics. 

Dr. H. Lowery, author of “The background of music’’, 
spoke of the difficulty of teaching music students who had 
only a negligible knowledge of science or mathematics, and 
said that he started with orchestration as the basis of his 
instruction. Lieutenant ROBERTS, who had been concerned 
with the teaching of young boys training to be army mu- 
sicians, also found the purely musical approach the most 
successful method. Dr. W. H. GEORGE (Chelsea Polytechnic) 
pointed out that musical acoustics was a very specialised 
subject; he thought that the theory of vibrations should be 
presented as a basic part of all physical theory and thereby 
do much to overcome the prejudice that a knowledge of 
music is essential to a study of acoustics. 

Professor FLoyp spoke of the necessity of giving pre- 
medical students a good understanding of the mechanism 
of the human ear, Dr. E. J. [Rons urged the use of graphs 
to illustrate wave-motion and Mr. S. H. Humsy cited ex- 
amples of how the ordinary architect was often years behind 
in his knowledge of the acoustic design of rooms. Other 
points emerging from later speakers were the correlation of 
acoustic with electromagnetic theory, and the possibility of 
making more space for acoustics in the school syllabus by 
incorporating the theory of harmonic motion and wave- 
motion in the applied mathematics course. 

A brief summary of the acoustics questionnaire sent to 
teaching institutions and industrial organisations by the 
Acoustics Group of the Physical Society was also presented 
by Dr. SrepHeEns. On the educational side it revealed that 
instruction in acoustics was a normal part of the B. Se. 
honours degree course in physics, but a few courses were 
also associated with equivalent degree examinations in 
music, architecture and engineering. Nearly 50 firms replied 
to the industrial questionnaire but only one or two salient 
features will be quoted here. It was surprising for instance 
to find only one out of seven firms in the acoustics material 
industry recruiting staff who have already received educ- 
ation in acoustics of a university degree standard. The 
electrical engineering and radio and telecommunication 
industries appear to be the chief industrial users of acoustic 
specialists, but it must not be overlooked that many find 
employment in the Post Office research section, the B. B.C. 
and Government research establishments. In the case of 
those industries who have research associations, it is the 
general practice for the various firms to pass over their 
acoustical problems to these associations. R.W.B.S. 


